NAVAL 

POSTGRADUATE 

SCHOOL 


MONTEREY,  CALIFORNIA 


DISSERTATION 


ROBUST  ADAPTIVE  CONTROL  USING  A  FILTERING 

ACTION 

by 

Jia-Horng  Yang 

September  2009 

Dissertation  Supervisor: 

Roberto  Cristi 

Approved  for  public  release;  distribution  is  unlimited 


REPORT  DOCUMENTATION  PAGE  |  Form  Approved  0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  respons  e,  including  the  time  for 
reviewing  instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and 
reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Washington  headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202^302,  and  to  the  Office  of 
Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503. 

I.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

_ September  2009 _ Dissertation _ 

4.  TITLE  AND  SUBTITLE:  Robust  Adaptive  Control  Using  a  Filtering  5.  EUNDING  NUMBERS 

Action _ 

6.  AUTHOR(S)  Jia-Homg  Yang  _ 

7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PEREORMING 

Naval  Postgraduate  School  ORGANIZATION  REPORT 

Monterey,  CA  93943-5000 _ NUMBER _ 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSORING  /  MONITORING 
N/A  AGENCY  REPORT  NUMBER 

II.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 

1 2a.  DISTRIBUTION  /  AVAII.ABII TTY  STATEMENT  I2b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited _ 

13.  ABSTRACT  (maximum  200  words) 

This  dissertation  describes  the  design  of  an  adaptive  controller  for  single-input  single-output  (SISO)  systems  with 
guaranteed  bounds  on  the  transient  response,  and  robustness  with  external  disturbances  and  unmodeled  dynamics. 
Developed  from  a  current  approach  called  “LI  adaptive  controller”,  we  show  that  by  adding  two  properly  designed 
low  pass  fdters  at  the  input  and  at  the  estimator  we  can  control  the  transient  response  and  the  sensitivity  of  the 
overall  system  to  external  disturbances  and  unmodeled  dynamics.  Global  stability  of  the  overall  adaptive  system  is 
mathematically  proven  under  the  assumption  that  the  system  is  minimum  phase  (i.e.,  with  the  zeros  of  the  transfer 
function  in  the  stable  region)  and  bounds  of  the  systems  parameters  are  known  to  the  designer. 

The  extension  of  this  approach  to  non-minimum  phase  systems,  such  as  systems  with  flexible 
appendages,  is  also  considered.  We  show  that  a  non-minimum  phase  plant  augmented  with  a  properly  designed 
parallel  system  results  in  a  minimum  phase  system.  The  augmenting  system  most  easily  comes  from  the  inverse  of 
a  stabilizing  Proportional-Integral-Derivative  (PID)  controller,  designed  to  be  least  sensitive  to  parameter 
uncertainties.  This  approach  is  applied  to  a  flexible  arm  in  a  testbed  at  the  Naval  Postgraduate  School,  called  the 
Flexible  Spacecraft  Simulator  (FSS),  which  emulates  realistic  conditions  in  space.  Experimental  results  prove  the 
effectiveness  of  the  controller  presented  in  this  dissertation. 


14.  SUBJECT  TERMS  15.  NUMBER  OE 

low  pass  filter,  LI  adaptive  controller,  unmodeled  dynamics,  non-minimum  phase,  PID  feedback,  PAGES 

flexible  problems  _ 121 _ 

16.  PRICE  CODE 

17.  SECURITY  18.  SECURITY  19.  SECURITY  20.  LIMITATION 

CLASSIEICATION  OE  CLASSIEICATION  OE  THIS  CLASSIEICATION  OE  OE  ABSTRACT 

REPORT  PAGE  ABSTRACT 

Unclassified  Unclassified  _ Unclassified _  UU 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 


1 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited 


ROBUST  ADAPTIVE  CONTROL  USING  A  FILTERING  ACTION 

Jia-Horng  Yang 
Major,  Taiwan  Army 

M.S.,  National  Defense  University  Chung  Cheng  Institute  of  Technology,  2000 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

DOCTOR  OF  PHILOSOPHY  IN  ELECTRICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September2009 


Author:  _ 

Jia-Horng  Yang 


Approved  by: 


Roberto  Cristi 

Professor  of  Electrical  and  Computer  Engineering, 
Dissertation  Supervisor  and  Committee  Chair 


Murali  Tummala  Xiaoping  Yun 

Professor  of  Electrical  and  Professor  of  Electrical  and 

Computer  Engineering,  Computer  Engineering, 


Gary  Hutchins  Donald  Wadsworth 

Associate  Professor  of  Electrical  and  Professor  of  Space  Systems 
Computer  Engineering,  Academic  Group, 

Approved  by:  _ 

Jeffrey  B.  Knorr,  Chair,  Department  of  Electrical  and  Computer  Engineering 

Approved  by:  _ 

Douglas  Moses,  Vice-Provost  for  Academic  Affairs 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


IV 


ABSTRACT 


This  dissertation  describes  the  design  of  an  adaptive  controller  for  single-input 
single-output  (SISO)  systems  with  guaranteed  bounds  on  the  transient  response,  and 
robustness  with  external  disturbances  and  unmodeled  dynamics.  Developed  from  a 
current  approach  called  “LI  adaptive  controller”,  we  show  that  by  adding  two  properly 
designed  low  pass  filters  at  the  input  and  at  the  estimator  we  can  control  the  transient 
response  and  the  sensitivity  of  the  overall  system  to  external  disturbances  and 
unmodeled  dynamics.  Global  stability  of  the  overall  adaptive  system  is  mathematically 
proven  under  the  assumption  that  the  system  is  minimum  phase  (i.e.,  with  the  zeros  of 
the  transfer  function  in  the  stable  region)  and  bounds  of  the  systems  parameters  are 
known  to  the  designer. 

The  extension  of  this  approach  to  non-minimum  phase  systems,  such  as 
systems  with  flexible  appendages,  is  also  considered.  We  show  that  a  non-minimum 
phase  plant  augmented  with  a  properly  designed  parallel  system  results  in  a  minimum 
phase  system.  The  augmenting  system  most  easily  comes  from  the  inverse  of  a 
stabilizing  Proportional-Integral-Derivative  (PID)  controller,  designed  to  be  least 
sensitive  to  parameter  uncertainties.  This  approach  is  applied  to  a  flexible  arm  in  a 
testbed  at  the  Naval  Postgraduate  School,  called  the  Flexible  Spacecraft  Simulator 
(FSS),  which  emulates  realistic  conditions  in  space.  Experimental  results  prove  the 
effectiveness  of  the  controller  presented  in  this  dissertation. 
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EXECUTIVE  SUMMARY 


Adaptive  control  systems  have  been  proposed  to  control  systems  with  uncertain 
dynamics,  widely  varying  under  different  operating  conditions.  They  have  been  the 
subject  of  research  for  a  number  of  years  and  stability  has  been  established  under  ideal 
conditions.  In  spite  of  the  amount  of  research  performed  on  this  class  of  control  systems, 
sensitivity  to  external  disturbances  and  modeling  errors  together  with  poor  transient 
response  performance  made  engineers  skeptical  to  their  applications.  In  particular,  the 
Model  Reference  Adaptive  Control  (MRAC)  approach,  which  seems  to  be  the  most 
effective  under  ideal  conditions,  can  easily  become  unstable  in  the  presence  of 
unmodeled  dynamics  not  accounted  by  the  system  model.  Where  instability  means  a 
complete  loss  of  the  platform,  such  as  in  space  applications,  adaptive  control  is  far  from 
being  considered  a  possible  solution  to  control  problems. 

In  this  dissertation,  we  address  some  of  these  issues,  such  as  transient  response 
and  robustness  in  the  presence  of  external  disturbances  and  modeling  errors.  The 
proposed  adaptive  control  system  is  derived  from  a  recently  introduced  adaptive  control 
approach  called  LI  Adaptive  Control,  with  LI  referring  to  the  norm  of  a  certain  operator 
affecting  the  tracking  error.  We  show  that  the  addition  of  two  simple  low  pass  filters, 
one  at  the  control  input  and  one  at  the  parameter  estimator,  eliminates  high  frequency 
components,  which  adversely  affect  both  the  transient  response  and  the  sensitivity  to 
disturbances.  These  two  additional  filters  must  be  properly  designed  in  order  to 
guarantee  stability  and  fast  convergence. 

A  major  drawback  of  this  proposed  scheme,  relative  to  more  traditional  adaptive 
systems,  is  the  need  for  a  nominal  model  for  the  system  and  bounds  on  the  parameter 
uncertainties.  In  other  words,  while  transitional  adaptive  controllers  assume  the  system  is 
almost  completely  unknown  (a  black  box),  in  the  proposed  scheme  we  assume  the 
system  to  be  within  bounds  (a  gray  box).  In  most,  if  not  all  applications,  the  gray  box 
assumption  is  more  realistic,  since  the  order  of  magnitude  of  the  dynamics  is  understood. 

A  further  issue  examined  is  the  extension  to  non-minimum  phase  systems.  This 
problem  is  addressed  by  replacing  the  actual  system  with  an  “augmented”  system, 

xiii 


implemented  by  adding  a  dynamie  system  in  parallel.  It  is  easy  to  show  that  a  non¬ 
proper  stabilizing  eompensator,  like  a  proportional-integral-derivative  (PID)  eontroller, 
ean  be  used  by  defining  the  inverse  of  its  transfer  funetion  (thus  proper)  as  the 
augmenting  system.  The  eombination  then  beeomes  minimum  phase,  and  it  ean  be 
eontrolled  by  an  adaptive  eontroller.  Applieation  to  a  flexible  arm  shows  the 
effeetiveness  of  this  result,  and  testing  at  the  Spaee  Researeh  and  Design  Center  (SRDC) 
at  the  Naval  Postgraduate  Sehool  emulates  a  flexible  roboties  arm  freely  moving  in  spaee. 

Although  all  stability  results  are  proven  mathematieally,  a  number  of  issues 
remain  to  be  addressed.  Speeifieally,  further  researeh  is  required  to  ehoose  the  optimal 
number  of  design  parameters  for  best  performanee  of  the  proposed  system. 
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I.  INTRODUCTION 


A.  MOTIVATION  AND  BACKGROUND 

The  particular  feature,  which  distinguishes  adaptive  control  systems  from  other 
control  techniques,  is  the  adaptive  loop  that  provides  real  time  estimates  of  the  Plant 
Dynamics. 

A  general  structure  of  an  adaptive  controller  is  shown  in  Figure  1 .  The  system  we 
want  to  control,  called  the  plant,  is  modeled  in  terms  of  a  set  of  parameters  9* ,  unknown  or 
partially  known.  The  goal  of  the  controller  is  to  estimate  the  plant’s  parameters  and  adjust 
the  controller’s  parameters  accordingly. 


el 


Figure  1.1  Model  Reference  Adaptive  Control. 
After  [1] 


1 


As  we  can  see  in  Figure  1,  there  are  two  loops.  One  is  the  standard  feedback  loop 
from  the  output  y  back  to  the  input  u  .  The  other  is  the  loop  from  the  input-output  of  the 
plant  to  the  compensator  parameters.  These  two  loops  can  interact  with  each  other  and 
cause  instability  in  the  system. 

This  approach  has  been  around  for  a  number  of  decades,  and  stability  under  ideal 
conditions  is  well  understood  [1]. 

In  spite  of  the  seemingly  attractive  feature  of  self-adjusting  to  time  varying 
dynamics,  researchers  and  engineers  have  always  been  skeptical  of  adaptive  control  and,  as 
of  today,  it  is  not  widely  used.  There  are  a  number  of  reasons  for  this  skepticism.  In 
particular: 

1.  Because  of  internal  modulation  effects,  fast  adaptation  produces  a  high 
frequency  input.  This  input  excites  unmodeled  dynamics  easily,  and  causes 
systematic  instability.  Typical  examples,  from  Rohrs  [2],  have  shown  that 
under  very  mild  conditions  of  unmodeled  dynamics  that  adaptive  control 
system  can  easily  become  unstable. 

2.  There  is  no  guarantee  on  the  bounds  of  the  transient  response  of  both  input  and 
output  signals.  Literature  review  [2-80]  shows  that  most  results  do  not  provide 
any  guidelines  on  transient  performance  analysis  or  explain  how  to  improve 
transient  performance. 

During  the  past  few  years,  there  has  been  greater  interest  in  adaptive  control,  as 
evidenced  by  its  prevalence  in  the  literature.  A  particular  approach  by  Naira 
Hovakimyan  [1]  shows  that  the  addition  of  a  filtering  term  in  the  loop  greatly  improves 
both  robustness  and  transient  response  behavior.  This  approach  is  called  “LI  adaptive 
controller”,  since  it  guarantees  an  arbitrarily  small  bound  on  the  transient  response 
using  the  small  gain  theorem  on  the  Lj  norms  of  the  error  operators  [81].  It  was 
developed  for  state  feedback  and  then  extended  to  input-output  systems.  However,  this 
extension  does  not  follow  the  standard  pole  placement  problem  in  classical  adaptive 
approaches  [82],  and  therefore,  yields  in  a  more  general  parameterization.  This  results 
in  an  adaptive  system  for  a  restrictive  class  of  linear  systems,  but  with  the  capability  of 
compensating  for  a  wider  class  of  nonlinearities. 

The  results  shown  in  [81,  83-88]  have  proved  the  effectiveness  of  this  approach 
in  terms  of  robustness  and  transient  behavior. 
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In  this  dissertation,  we  address  the  problems  of  designing  a  elass  of  Adaptive 
Control  systems  whieh  yield  fast  adaptation,  thus  good  transient  response,  and 
robustness  to  modeling  errors. 

In  addition,  we  show  that  this  approaeh  ean  be  extended  to  non-minimum  phase 
systems,  where  the  eurrent  approaehes  are  unable  to  stabilize  the  system.  Although  this 
approaeh  requires  more  knowledge  about  the  system  in  order  to  eontrol  it,  it  is  still 
attraetive  in  oases  where  preoise  eontrol  is  needed  in  the  presenoe  of  model 
unoertainties.  The  typioal  ease  we  analyze  is  the  preoision  eontrol  of  a  flexible  arm  in 
spaoe. 

B,  RESEARCH  GOALS 

After  introduoing  the  properties  of  a  general  adaptive  eontrol  system  and  the 
Model  Referenoe  Adaptive  Control  (MRAC)  system,  we  examine  the  drawbaoks  of 
these  systems  and  use  LI  adaptive  eontrol  to  address  them.  With  the  proposed  LI 
adaptive  oontroller,  not  only  is  the  traoking  error  driven  to  zero,  but  the  transient 
response  is  eontrolled,  due  to  adaptation.  The  final  result  is  a  eontroller  whieh,  under 
the  assumption  of  low  disturbanees,  yields  fast  adaptation  to  parameter  variations. 

Another  goal  is  to  strengthen  the  robust  ability  of  the  adaptive  system.  After 
proper  ehoiee  of  low  pass  filter,  referenee  model,  and  other  important  elements,  the 
proposed  adaptive  eontroller  demonstrates  an  improvement  in  the  system’s  ability  to 
deal  with  the  problems  of  modeling  errors,  nonlinearities  and  external  bounded 
disturbanees. 

C.  DISSERTATION  ORGANIZATION 

This  dissertation  is  organized  as  follows:  Chapter  II  introduees  the  fundamentals  of 
adaptive  eontrol  systems,  speeifieally  the  direet  Model  Referenee  Adaptive  Control 
(MRAC)  approaeh.  Rather  than  repeating  the  detailed  mathematieal  derivations,  whieh 
are  found  in  the  standard  referenees,  a  simple  first  order  example  is  presented,  whieh  has 
all  the  “fiavor“  of  a  full  theoretieal  presentation,  but  does  not  inelude  the  “fat”.  Also  in 
this  ehapter,  analytieal  arguments  and  eomputer  simulations  are  used  to  demonstrate  that 
fast  adaptation  results  in  high  frequeney  signals  and  a  possible  loss  of  stability. 
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Additionally,  using  this  same  example,  it  is  shown  that  the  addition  of  a  simple  low  pass 
filter  in  the  loops  guarantees  good  transient  performanee  and  better  robustness  to 
modeling  errors. 

Chapter  III  diseusses  the  general  extension  to  a  system  of  arbitrary  degree. 
Typieally,  the  plant  in  a  MRAC  system  must  be  minimum  phase.  This  means  that  the 
plant  to  be  eontrolled  ean  be  unstable,  and  that  the  zeros  of  its  transfer  funetions  must  be 
in  the  left  half  side  s-plane  (in  the  “stability  region”).  Additionally,  the  use  of  a  priori 
knowledge  of  plant  nominal  dynamies  and  its  uneertainties  is  diseussed  in  the  design  of 
an  adaptive  eontroller,  which  is  globally  stable  with  a  transient  response  having  an 
arbitrarily  small  bound.  This  response  is  obtained  provided  the  effect  of  external 
disturbances  is  negligible. 

Chapter  IV  addresses  the  minimum  phase  assumption,  where  the  adaptive  system 
is  extended  to  a  class  of  non-minimum  phase  systems.  It  is  shown  that  a  non-minimum 
phase  system  can  be  “augmented”  to  a  minimum  phase  transfer  function  by  adding  a 
properly  designed  parallel  system  whose  response  is  zero  at  steady  state.  Further,  it  is 
shown  that  if  the  system  is  stabilizable  by  a  standard  proportional — integral — derivative 
(PID)  controller,  then  the  transfer  function  of  the  added  system  is  just  the  inverse  of  the 
PID  controller.  This  is  particularly  attractive,  since  a  properly  design  PID  controller  is 
generally  robust  in  the  presence  of  model  uncertainties.  However,  such  a  system  has  the 
drawback  that  its  transfer  function  is  not  proper,  due  to  the  derivative  action  of  the 
controller.  But,  since  its  inverse  is  implemented  in  conjunction  with  the  adaptive 
controller,  its  transfer  function  becomes  strictly  proper,  and  thus  easily  realizable.  This 
chapter  also  discusses  an  experimental  implementation  of  this  approach  on  a  flexible  arm 
in  the  satellite  Research  and  Design  Center  (SRDC)  at  the  Naval  Postgraduate  School 
(NPS).  In  this  experiment,  vibrations  due  to  highly  damped  flexible  modes  in  the  system 
are  properly  compensated  for  so  that  the  position  of  the  end  effector  quickly  settles  to  the 
desired  position. 
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Finally,  conclusions  and  further  research  are  presented  in  the  final  ehapter  of  this 
dissertation.  Additionally,  in  order  to  better  illustrate  the  proposed  adaptive  algorithm,  a 
number  of  computer  simulations  and  the  results  of  applieations  are  shown  alongside  the 
theoretieal  presentation. 


D,  NOTATION 

In  this  dissertation,  we  define  “  p  ”  the  differential  operator  as 

px{t)  =  x{t) 


Also,  the  notation 


y{t)  = 


B{p) 

A(p) 


u(t) 


with  B(p) ,  A(p)  polynomials  in  p  ,  indieate  the  differential  equation 

p^y{t)  +  a^p^^^yit)  +  •  •  •  +  %y(0  =  b^p^^^u{t)  +  •  •  •  +  b^M{t) 
with  a. ,  bj  the  eoeffieients  of  A{p) ,  B{p) ,  respeetively. 

In  all  situations  where  the  Laplaee  Transform  exists,  the  differential  operator  p  ean  be 
also  interpreted  as  the  Laplaee  variable  s. 
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II.  INTRODUCTION  TO  ADAPTIVE  CONTROL 


A.  PREFACE 

According  to  Webster’s  dietionary,  to  adapt  means  “to  ehange  (oneself),  so  that 
one’s  behavior  will  eonform  to  new  or  ehanged  eireumstanees.”  This  is  exactly  the  goal 
of  an  adaptive  eontrol  system,  to  self  adjust  to  ehanging  dynamies  and  environmental 
eonditions.  A  typieal  example  is  the  design  of  a  high-performanee  jet  fighter  or  a  fire 
helieopter  [1],  They  both  operate  over  a  wide  range  of  speeds,  altitudes,  and  weights 
(before  and  after  loading  water  for  the  ease  of  a  lire  helieopter)  and  their  dynamies  are 
nonlinear  and  time  varying  as  deseribed  by 

x{t)  =  A{t)x{t)  +  B{t)u{t) 
y{t)  =  C{t)x{t)  + D{t)u{t) 

where  x(0)  =  x^  is  assumed  to  be  the  initial  state,  and  A{t)  ,  B{t) ,  C{t) ,  and  D{t)  are 
funetions  of  the  operating  time  t. 

As  the  system  goes  through  different  operating  eonditions,  the  dynamie  model 
A{t)  ,  B{t)  ,  C{t)  ,  and  D{t)  ehanges  eonsiderably.  The  adaptive  eontroller  strueture 
eonsists  of  a  feedbaek  loop  and  a  eontroller  with  adjustable  gains  as  shown  in  Figure  2.1. 


Figure  2.1.  Controller  Strueture  with  Adjustable  Controller  Gains. 

From  [1] 


7 


In  this  chapter,  two  different  approaehes  to  adaptive  eontrol  systems  are  briefly 
introdueed.  They  are  the  Indireet  Model  Referenee  Adaptive  Control  (Indireet  MRAC) 
and  the  Direet  Model  Referenee  Adaptive  Control  (Direet  MRAC)  approaehes.  This  is 
followed  by  a  simple  first  order  example  to  show  the  teehniques  to  analyze  eonvergenee 
and  stability.  Finally,  the  approaeh  presented  in  this  researeh  is  introdueed  for  the  first 
order  example  to  demonstrate  the  main  features  of  the  proposed  algorithm. 

B,  MRC,  INDIRECT  MRAC,  AND  DIRECT  MRAC 

MRAC  is  derived  from  the  Model  Referenee  Control  (MRC)  approaeh.  A  linear 
time  invariant  (LTI),  SISO  plant  of  MRC  is  shown  in  Figure  2.2. 


Figure  2.2.  Model  Referenee  Control. 

From  [1] 

In  MRC,  full  knowledge  of  the  plant  dynamies  is  assumed.  The  eontroller  must 
be  designed  to  meet  the  designer’s  requirements,  identified  in  the  referenee  model,  as  the 
desired  Input-Output  (I/O)  properties  of  the  elosed  loop.  The  goal  of  the  MRC  approaeh 
is  to  find  a  proper  feedbaek  eontrol  law  so  that  the  elosed  loop  matehes  a  referenee  model 
with  transfer  funetion  Wm(p). 

The  feedbaek  eontroller  C{0*^  is  designed  so  that  all  signals  are  bounded  and 
the  elosed  loop  transfer  funetion  from  r{t)  to  y{t)  is  equal  to  W^{p) .  This  implies  that 
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the  object  of  the  desired  controller  is  to  drive  the  tracking  error 

(0  =  k(0  “  km  (0  close  to  zero.  The  requirements  of  the  MRC  are  as  follows: 

1 .  The  zeros  of  the  plant  must  be  minimum  phase. 

2.  A  good  understanding  of  the  plant  and  performance  is  required. 

If  the  zeros  of  the  plant  are  non-minimum  phase,  the  MRC  may  have  an  unstable  pole 
zero  cancellation  thus  causing  instability.  When  the  plant  parameters  are  not  totally 
known,  an  adaptive  system  based  on  estimates  of  the  plant  parameters  must  be  developed. 
The  estimated  parameters  can  be  obtained  by  direct  and  indirect  approaches.  The 
resulting  control  schemes  are  known  as  MRAC.  A  block  diagram  of  a  direct  MRAC 
system  is  shown  in  Figure  2.3,  and  an  indirect  MRAC  system  is  shown  in  Figure  2.4. 

The  difference  between  a  direct  MRAC  and  indirect  MRAC  is  that  a  direct 
MRAC  updates  the  estimated  controller  parameter  vector  directly  from  input  and  output 
data.  In  indirect  adaptive  control,  the  controller  parameters  are  computed  on  the  basis  of 
the  estimated  plant  dynamics.  In  this  dissertation,  the  approach  is  based  on  the  direct 
MRAC.  Accordingly,  the  concept  of  the  direct  MRAC  is  discussed  in  the  next  section. 
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F igure  2.3.  Indirect  MRAC . 

After  [1] 


Input 

Comma] 

r{t) 


Figure  2.4.  Direct  MRAC. 
From  [I] 
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c. 


MODEL  REFERENCE  ADAPTIVE  CONTROL 


1,  Introduction  of  Model  Reference  Adaptive  Control 

The  general  idea  behind  MRAC  is  to  ereate  a  elosed  loop  eontroller  with 
estimated  parameters  that  ean  be  updated  direetly  from  input/output  observations.  A 
standard,  simple  direet  MRAC  is  shown  in  Figure  2.5. 


As  shown  in  Figure  2.5,  the  output  of  the  system  y{t)  is  eompared  to  a  desired 
response  y„{t)  from  a  referenee  model.  The  eontrol  parameters  are  updated  based  on 
the  error  e{t)  between  y{t)  and  y^{t)  .  The  goal  of  the  MRAC  is  to  update  the 
estimated  parameters  d{t)  so  that  the  plant  response  matehes  the  response  of  the 
reference  model  (p)  .  In  what  follows,  a  first  order  example  of  an  adaptive  control 
system  is  presented  and  the  methodology  for  determining  stability  and  convergence  is 
shown. 


11 


2,  Theory  of  Model  Reference  Adaptive  Control 

A  simple  example  of  a  first  order  case  is  used  to  introduce  the  theory  of  MRAC. 
Assume  a  plant  with  transfer  function 

W^{p)  =  -^  (2.1) 

p-a 

where  a  is  unknown,  and  it  is  not  necessarily  positive  or  negative.  This  implies  that  no 
assumption  regarding  the  original  plant’s  stability  is  made.  Equation  (2.2)  depicts  the 
state  space  model  as 


x{t)  =  ax{t)  +  u{t) 
y{t)  =  x{t). 

Let  the  reference  model  have  dynamics 

(0  =  +  v(t) 

where  a^>0  for  stability.  The  reference  output  signal  (t)  represents  the  desired 
output  the  plant  has  to  track.  The  goal  of  the  controller  in  MRAC  is  to  drive  the  tracking 
error  e{t)  =  y{t)-y^{t)  close  to  zero.  In  order  to  do  so,  the  plant  dynamics  represented 
by  Equation  (2.2)  are  expressed  in  terms  of  the  reference  model  as 


x(t)  =  -a^^x(t)  +  (a„  +  a)x(t)  +  u(t).  (2.4) 

Eurther,  the  unknown  parameter  is  defined  as 


6  =  a^^+a  (2.5) 

so  that  Equation  (2.4)  can  be  written  as 

x{t)  =  -a^x{t)  +  6x{t)  +  u{t).  (2.6) 

This  yields  the  structure  of  the  control  input,  in  terms  of  the  model  parameter  6 ,  as 
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u{t)  =  -6x{t)  +  v{t) 


(2.7) 


which  yields  the  desired  closed  dynamies 

x{t)  =  -a^x{t)  +  v{t).  (2.8) 

If  we  know  the  parameter  a  in  equation  (2.2),  then  we  can  calculate  the  control 
parameter  in  equation  (2.5),  and  use  the  control  input  signal  as  in  equation  (2.7)  to  obtain 
the  desired  eontrol  behavior.  However,  since  we  do  not  know  the  system  dynamics,  we 
use  a  time  varying  estimated  parameter  6{t)  is  used  to  replace  the  true  eontrol  parameter 
6 .  In  this  way,  the  eontrol  becomes 

u{t)  =  -6{t)x{t)  +  v{t).  (2.9) 

Now,  the  problem  is  how  to  compute  the  estimated  parameter  6{t)  based  on  input/output 
data.  In  order  to  do  this,  the  eontrol  input  expression  of  Equation  (2.9)  is  substituted  into 
the  plant  model  Equation  (2.6)  resulting  in. 

x{t)  =  -a^^x{t)->r6{t)x{t)  +  v{t)  (2.10) 

where  the  parameter  error  6{t)  =  6-  6{t)  is  the  difference  between  actual  value  6  and 
estimated  parameter  6{t)  .  Now  if  the  tracking  error  e{t)  is  defined  as 

e{t)  =  y{t)  -  (t)  =  x{t)  -xjt)  (2.11) 

and  Equations  (2.10)  and  (2.3)  are  eombined  to  relate  the  parameter  error  0{t)  with  the 
tracking  error  e{t) ,  the  result  is 

e{t)  =  -a^^e[t)  +  d{t)x{t).  (2.12) 

In  order  to  provide  a  viable  estimate,  the  parameter  0{t)  is  updated  as  follows 
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(2.13) 


6{t)  =  -6{t)  =  /x{t)e{t) 

where  /  >  0  is  an  arbitrary  eonstant.  The  basis  for  this  is  shown  as  follows. 

Define  the  Lyapunov  funetion 

K(e(0,^(0)  =  i|e(0f  +^r''  \e{t)\  (2.14) 

Differentiating  V  along  the  trajeetories  of  the  systems  (2.12),  (2.13),  and  (2.14),  we 
obtain 

V (e(0, 0{t))  =  e{t)i^-a^e[t)  +  6{t)x{t))  +  / '^(t) (-/^(OKO)  (2.15) 

whieh  yields 

V{e{t)M))  =  -aM^t  <0  (2.16) 

By  standard  argument  this  shows  that  the  estimated  parameter  is  bounded  and  the 
error  goes  to  zero  [1]. 

In  this  example,  we  ean  also  verify  that  an  upper  bound  of  the  error  e{t)  ean  be 
easily  determined.  Infaet,  sinee  V{e{t),6{t))  is  non  inereasing 

V{e{0),m)>V{e{t)M))  (2.17) 

we  ean  write  from  (2.14)  that 

i|e(0)|"  +  iy-'  |#(0)f  >  |e(Or  +'-r'  |^(of  ■  (2.18) 

By  ehoosing  zero  initial  eondition  x^(0)  =  x(0)  ,  i.e.,  e(0)  =  0  we  obtain 

i r'  |0(O)f  >  t |e(0f  +^r-'  \kt)\  >  i|e(f)f  (2.19) 

for  all  t  >  0  .  This  shows  an  upper  bound  for  the  traeking  error  as 

,  ,2  ^(0)' 

\e{t)\  < -  (2.20) 

y 
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for  all  ^  >  0  .  From  Equation  (2.20),  the  relationship  between  tracking  error  e{t)  , 
adaptive  gain  y ,  and  initial  parameter  error  ^(0)  is  found.  By  increasing  the  gain  y ,  the 
maximum  error  can  be  made  arbitrarily  small,  thus  achieving  any  desired  tracking. 
Figures  2.6  and  2.7  show  simulations  for  different  values  of  the  adaptation  gain  y . 


input  u(t) 

2  p  ■  ^  ^  I  I  I  I  I 

0  - 

-2  ^ 

_4  I _ I _ 1 _ 1 _ 1 _ I _ I _ 1 _ 

0  5  10  15  20  25  30  35  40 


sec 

plant  x(t)  and  model  x^^it)  outputs 


Figure  2.6.  Response  of  Modeled  MRAC  Adaptive  System  with  y  =  50  . 
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Figure  2.7.  Response  of  Modeled  MRAC  Adaptive  System  with  y  =  1000 . 


These  simulations  show  that,  by  inereasing  the  adaptation  gain  y ,  the  maximum 
traeking  error  deereases.  Flowever,  this  eomes  with  a  price.  As  shown  in  the  simulation, 
the  input  signal  exhibits  a  higher  frequency  component.  In  the  next  section,  these  high 
frequency  components  will  be  discussed,  including  their  excitation  of  unmodeled 
dynamics,  which  cause  the  system  to  lose  tracking  or  to  become  unstable.  This  is  seen  in 
Figures  2.8  and  2.9,  where  the  addition  of  unmodeled  poles  at  /»  =  -20  +  7*5  and 

425 

V  =  -20  -  /  *  5  ,  with  transfer  function  — r -  ,  causes  the  system  to  either 

/ +40^7  +  425 

becomes  unstable  or  loose  tracking  (Figure  2.9). 
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Figure  2.8.  Response  of  Modeled  MRAC  Adaptive  System  with  y  =  50  . 
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Figure  2.9.  Response  of  MRAC  with  Unmodeled  Dynamics  and  y  =  50  . 
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Figure  2.9  shows  that,  by  increasing  the  adaptation  gain  y  in  the  MRAC,  the 
maximum  tracking  error  may  not  decrease  due  to  the  unmodeled  dynamics  thus  leading 
to  loss  of  tracking. 

3.  Adaptive  Control  with  a  Filtering  Action 

The  arguments  presented  in  Chapter  II.  C.2  show  that  there  between  faster 
adaptation  (large  y)  comes  at  the  expenses  of  a  high  frequency  transient  response 
which  can  cause  instabilities  in  the  presence  of  modeling  errors.  This  can  be  seen 
from  the  input  signal,  repeated  here  for  convenience  as, 

u{t)  =  -6{t)x{t)  +  v{t). 

Using  the  definition  of  parameter  error  6{t)  =  6 -6{t)  and  its  relation  to  the  tracking 
error  e{t)  shown  in  Equation  (2.12),  the  input  signal  can  be  expressed  as 

u{t)  =  -6x{t)  +  v{t)  + e{t)  + a^e{t).  (2.21) 

While  the  error  e{t)  can  be  made  arbitrarily  small  by  increasing  the  adaptation  gain  y , 
its  derivative  e{t)  not  only  is  not  necessarily  small,  but  it  can  become  quite  large  due  to 
the  high  frequencies  in  the  transient  response. 

A  better  overall  response  is  obtained  by  introducing  a  Low  Pass  Filter  with 
transfer  function  (5)  at  the  input  signal,  which  yields 

u{t)  =  C,{p)[-e{t)x{t)  +  v{t))  (2.22) 

This  filter’s  bandwidth  should  be  larger  than  the  reference  model’s,  so  that  the  filter  stops 
the  high  frequencies  generated  by  the  adaptive  system  without  affecting  the  dynamics  in 
steady  state. 

The  input  signal  u{t)  can  be  written  again  in  terms  of  6{t)  and  e{t)  as 

u{t)  ^  Cp {p)[-ex{t)  +  v(0)  +  Cp {p){e{t)  +  a„,e(0) •  (2.23) 
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The  rightmost  term  shows  the  effect  of  the  adaptation  error,  and  it  can  be  written  as 

u{t)  =  [p  +  a)Cp{p)e{t)  (2.24) 

where  the  transfer  function  (5  +  a)C^(5)  is  proper  (i.e.,  no  differentiation)  and  stable. 
This  implies  that  small  values  for  e{t)  also  result  in  small  values  for  u{t)  ,  and  the 
controller  is  always  close  to  an  ideal  controller  with  parameter  6 . 

Figures  2.10  and  2.11  show  this  behavior.  The  system  is  the  same  as  the  system 
presented  in  Chapter  II.  C.2,  but  it  is  observed  that  an  increased  adaptation  gain  y  yields 
faster  response,  with  better  transient  behavior.  This  adaptive  control  with  a  filtering 
action  approach  has  been  called  LI  adaptive  control  [81,  83-88]. 


Figure  2.10.  Response  of  LI  Unmodeled  MRAC  Adaptive  System  with  y  =  50 . 
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Figure  2.11. 


Response  of  LI  Unmodeled  MRAC  Adaptive  System  with  y  =  1000 . 
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III.  LI  ADAPTIVE  CONTROL  FOR  MINIMUM  PHASE  SYSTEMS 

A.  OVERVIEW  OF  LI  ADAPTIVE  CONTROLLER 

In  this  chapter,  the  adaptive  control  formulation  is  extended  to  a  general  SISO 
system.  In  particular,  it  is  shown  that  the  plant  ean  be  parameterized  in  terms  of  the 
eontroller  coeffieients,  so  that  they  ean  be  estimated  direetly  from  the  input  and  output  data. 
Further,  the  effeetiveness  of  adding  a  LPF  in  the  controller  is  addressed.  By  proper  ehoiee 
of  this  LPF,  it  is  shown  that  the  magnitude  of  the  transient  response  due  to  adaptation  ean 
be  made  arbitrarily  small  for  both  input  and  output  signals.  The  cost  is  a  perturbation  of  the 
reference  model,  which  is  still  globally  stable. 

This  chapter  is  organized  as  follows.  Model  definitions  are  introduced  in 
Chapter  III.  B.  Parameter  estimation  is  given  in  Chapter  III.  C.  Then,  all  properties  and 
elements  are  combined  to  show  the  overall  model  in  Chapter  III.  D.  In  Chapter  III.  E, 
the  LI  adaptive  controller  is  analyzed  under  ideal  conditions,  without  any  disturbances 
or  modeling  error.  Then,  in  Chapter  III.  F,  analysis  of  the  ideal  case  is  extended  to  the 
eonditions  of  external  disturbance  and  unmodeled  dynamies.  Finally,  in  Chapter  III.  G 
some  examples  are  eonsidered  that  eompare  MRAC  systems  and  verify  the  theorem  of 
LI  adaptive  controller.  After  this,  a  deeper  exploration  of  the  properties  is  possible. 

B,  MODEL  DEFINITIONS 


In  the  SISO  cases,  eonsider  a  system  with  dynamics  described  by 


yit)  =  K  w(0  +  Wo(0 

A{p) 

(3.1) 

where 

Ap)  =  P^+  + ...  + 

(3.2) 

B{p)  =  p^  +b,p^-A...  +  b^ 

(3.3) 

are  monie.  Then  the  referenee  model  is  defined  as 
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(3.4) 


where 


ymit)  = 


1 

D{p) 


v{t) 


D{p) = +...+d 


N-M 


(3.5) 


is  a  Hurwitz  polynomial  with  degree 

dD{p)  =  dA(p)  -  dB{p)  (3 .6) 

and  v{t)  is  an  arbitrary  external  bounded  input. 

In  the  aforementioned  diseussion,  and  in  most  of  the  related  literature,  the 
assumptions  on  the  plant  are  as  follows: 

a.  The  degree  of  the  numerator  B{p) ,  M,  and  the  denominator  A{p) ,  N,  are 
known. 

b.  The  sign  of  the  high  frequeney  eoeffieient  is  known  (assume  ^0  >  0 
without  loss  of  generality). 

e.  The  high  frequeney  eoeffieient  and  the  parameters  bf^, 

of  the  transfer  funetion  of  the  plant  are  unknown. 

d.  The  plant  ean  be  stable  or  unstable.  However,  the  numerator  B{p)  is 
Hurwitz,  i.e.,  it  has  roots  with  negative  real  parts. 

1.  Plant  Parameterization 

Now,  the  goal  is  to  parameterize  the  system  in  terms  of  the  eontroller  parameters. 
In  partieular,  following  the  pole  plaeement  approaeh,  and  assuming  the  following  state 
spaee  model,  as  follows: 

x(t)  =  Ax(t)  +  BkaU(t) 

(3.7) 

y{t)  =  C^x{t) 
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Using  an  observer  and  state  feedbaek,  the  closed  loop  poles  can  be  arbitrarily 
placed  as  shown  in 

x{t)  =  (A^  -  K^CJx(t)  +  B^kou(t)  +  K^y(t) 
kgU(t)  =  -LpX{t)  +  v{t). 

If  the  controller  gain  is  chosen  so  that 

det(5/  -A^+  B^L^)  =  B{p)D{p)  (3.9) 


and  the  observer  gain  as 


dQi{sI-A^+Kf^)  =  P{p) 


(3.10) 


with  P{p)  arbitrary,  then  the  Hurwitz  closed  loop  dynamics  becomes 


y{t)  = 


B{p) 


D{p)B{p) 


v{t) 


1 

D{p) 


v{t). 


(3.11) 


In  order  to  determine  a  suitable  plant  parameterization,  the  controller  in  Equation 
(3.8)  is  written  as 

k,u{t)  =  -L^  {si  -A^+  r  B^u  (0  -  {si  -A^+  )'  K^y{t)  +  v{t). 


This  expression  can  also  be  written  as 


k^u{t)  =  - m(0  -  y{t)  +  v{t) 
“  P{p)  P{p) 


where 


”  P  P  P'  P  p(^p)  p(^) 


P{p)  P{p) 


(3.12) 


(3.13) 

(3.14) 
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Define  the  veetors  of  eoeffieients 


and  the  filtered  signals  as 


f-]  P{p) 


P{p) 


N-2 

P 

P 

1 

jA'-l 

,A'-2 


P 

1 


u{t) 


y{t) 


then  the  plant  model  ean  be  written  as 


1 


yi^)  =  P^iKu{t)  +  h  (p,it)  +  k  ^(0). 
Dip)  -  ^ 


Sinee  the  bounds  on  k^  are  assumed  to  be  known  as 


equation  (3.19)  ean  be  expressed  as 


1 


yit)  =  +  [k^-k^)uit)  +  /z  (0  +  k''  <ty  (0) 

In  summary,  the  following  statement  ean  be  made: 

Given  the  plant 

yit)  =  K^jf\u{^) 

Aip) 


(3.15) 

(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 
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and  an  arbitrary  Hurwitz  polynomial  D{p)  with  degree  dA(p)  -  dB(p) ,  the  plant  ean  be 
represented  as 


=  r-j— 

Dip) 


(3.21) 


where 


(pit)  = 


u(t). 


N-\ 


Pip) 


uit), 


1 


N-\ 


-u{t) 


Pip)  Pip) 


yit),- 


Pip) 


yit) 


(3.22) 


and  6^  e  ^27v+i  controller  parameters. 

It  can  be  easily  seen  from  Equation  (3.21)  that,  in  order  to  track  the  model  — — , 

Dip) 

the  controller  has  to  be  of  the  form 

kjiit)  =  -^  (l)it)  +  vit)  (3.23) 


C.  PARAMETER  ESTIMATION 

In  this  section,  the  problem  of  estimating  the  parameters  ^  from  the  plant 
dynamics 

yit)  =  -  j—  ikjiit)  +  ^^^(0)  +  wit)  (3 .24) 

Dip) 

is  addressed.  In  order  to  do  so,  let  Cj^ip)  be  a  transfer  function  such  that  the  C^ip)Dip) 
is  proper  and  is  of  the  form 


C,iP)  = 


is  +  a)i,Y 


(3.25) 


with  L>  N -M  .  Thus  equation  (3.24)  can  be  arranged  as 

Ce  ip)Dip)yit)  =  Q  ip)ik„uit)  +  (kg  -  )uit)  +  h^(p^it)  +  k^  (j)^,  it))  +  it)  (3 .26) 


or 


CEip)Dip)yit)  =  C^ip)ikjiit)  +  ^(l)it))  +  w^it) 
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(3.27) 


where  w,,{t)  =  Cj,{p)D{p)w{t) .  If  w,;{t)  is  set  as  the  upper  bound  of  Wj,{t) ,  then  the 
parameter  of  equation  (3.26)  and  (3.27)  ean  be  written  as: 

CE{p)D{p)y{t)  -  {t)u{t)  =  e^X  (0  +  (t)  (3 .28) 

where  0^  e  and  X{t)  =  Cj,{p)(l){t)  . 

Notiee  that  the  parameter  veetor  0  appears  in  (3.28).  Therefore,  in  order  to 
estimate  it,  eall  ^(t)  an  estimate  and  define  the  error  as 

e{t)  =  C,{p)D{p)y{t)  -  Kp,{t)u{t)  -  6^  {t)X{t)  (3.29) 

This  ean  be  easily  related  to  the  parameter  error  6{t)  =  6- 6{t) 

e{t)  =  e'^{t)X{t)  +  w^{t).  (3.30) 


Based  on  this  equation,  the  parameter  6{t)  ean  be  estimated  as  follows: 

xxt) 


a  w^{t)  +  X{t) 


-  e{t),  /  =  0, . . . ,  IN 


(3.31) 


with 


0,  if  6».(0  =  0-^^  and  X.{t)e{t)  >  0 


Pi  (0  =  T  0,  if  (t)  =  6»,  and  X,.  {t)e{t)  <  0 


^  p ,  otherwise 


(3.32) 


and  p  >  0  arbitrary.  By  this  definition,  the  estimated  parameters  0-{t)  is  bounded  as 
shown  in  Figure  3.1.  When  the  estimated  reaehes  its  maximum  (minimum)  value, 
adaptation  is  turned  off  (//,.(t)  =  0)  as  long  as  it  tends  to  inerease  (deerease). 
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Figure  3.1.  The  Figure  of  Estimated  Parameter  0-  (t) . 


From  the  aforementioned  information,  it  can  be  shown  that  for  all  t>0  that 


r 


a-  o.^^<ext)<e,^ 


< 


b. 


0^  {t)X (0  <  s,  (0  ( cr  I (0|  +  ||X(0||) 


c. 


0{t) 


<lusXt) 


V. 


with  Sy  (t) ,  ^2  (0  >  such  that  for  all  t,  T  >  0 


(3.33) 
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^  "eaf 


t+T 


\\sXr)\dT<  < 


2/u 


0{t) 


if  Wj,  (r)  =  0  for  all  r 


(3.34) 


^ T  otherwise 


^  2//  o- 

The  proof  is  shown  in  Appendix  A. 

In  the  absence  of  disturbances,  (t)  =  0  for  all  t  >  0  .  Therefore,  Equation 
(3.31)  can  be  restated  as 


^(0  =  //,(0^^e(0,  /  =  0,...,2A 
||^(0|| 


(3.35) 


with 


^  0,  if  e,{t)  =  and  X.{t)e{t)  >  0 


//.(O  =  S  0,  if  e^t)  =  and  X.{t)e{t)  <  0 


(3.36) 


ju ,  otherwise 


In  the  case  of  no  disturbance,  in  order  to  avoid  singularities  during 
implementation,  the  denominator  ||A(t)||^  can  be  replaced  with  a  +  ||A(t)||^  for  any 
arbitrary  constant  a ,  without  loss  of  generality.  In  developing  the  simulated  model, 

II  l|2 

the  arbitrary  constant  is  set  to  one,  a  =  1 .  Thus,  the  denominator  becomes  1  +  ||A(t)|  . 

D,  OVERALL  MODEL 

As  in  the  first  order  system  example  presented  in  Chapter  II,  the  filter  {p)  is  a 
low  pass  filter  with  a  DC  value  of  one,  C^(0)  =  1 .  This  filter  is  critical  to  the  stability 

and  performance  of  the  overall  system,  and  it  is  addressed  in  greater  detail  later  in  this 
section,  including  the  criteria  for  its  selection. 

Let  the  input  of  the  adaptive  control  system  be  represented  as 
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(3.37) 


where  v(^)  is  the  external  input  and  ^(t)  is  as  previously  described  in  Equation  (3.22). 
The  estimated  parameters  0.{t)  are  computed,  as  shown  in  Equation  (3.31). 

Now  the  control  input  in  (3.37)  is  combined  with  the  parameterized  plant  model 
in  Equation  (3.24)  to  obtain  the  closed  loop  dynamics  model  as 

y(t)  =  v(t))  +  (e'm)  +  «(t) 

Dip)  D{p)  ^  ^  (3.38) 

k^u{t)  =  Cp  ip)[0^  it)(pit)  +  v(0)  -  Cp  ip)[0^(pit)) 

where  0^  =[^o“^m  the  true  parameters  vector.  Again  0{t)  ,  the  estimate, 

is  computed  as  in  Equation  (3.31),  and  0{t)  =  0- 0{t)  is  the  parameter  error. 


egit) 

Eigure  3.2.  Dynamics  Model  of  the  Adaptive  System. 


In  order  to  better  understand  the  global  stability  of  the  overall  adaptive  system, 
equation  (3.38)  can  be  represented  by  the  dynamic  model  in  Eigure  3.2  where  we  define 
the  error  term 
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eg{t)  =  Cp{p)(e^{t)^{t)) 


(3.39) 


as  the  effect  of  parameter  error.  The  overall  model  obtained  by  combining 

equation  (3.38)  with  the  two  plant  parameterizations 

A(p) 

and 


y(0  =  r-j—  + ^^^(0) 

D(p) 

where  .  This  yields  the  dynamics  model 


t(0 

M(t) 


Href  [p)Cf {p)[6\t)(l){t)  +  v(0)  +  H^[p) Wait) 


(3.40) 


defined  as  follows: 

y(t)  =  Hyy{p)y(t)  +  Hy^{p)Cp(p)[0^(t)<f>(t)  +  v(t))  +  H^^{p)wa(t) 

kjiit)  =  {p)y(t)  +  //„„  {p)Cp(p)[0^(t)<f>(t)  +  v(0)  +  //„„  {p)  Wait) 


where 
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H„(p)=(\-C,(p))\\- 


A{p) 


B{p)D{p) 


\  +  YCp[p) 


HyAp)  = 


(1+/) 


[\+rC,{p))D{p) 


^Ap) 


(l-C^(;?))  f  A(p)  ^  k(p) 

(\  +  rCj,{p))[B{p)D{p)  P{p)D{p) 


\ 

1  +1 

y 


k{p) 

P{p)D{p) 


HuAp) 

Huy{p) 

Hu«{p) 


1 

\  +  YCp{p) 

Cf{p)D{p)(  A{p)  ^ 

\  +  YCp{p)  yB{p)D{p)  j 

Cf{p)D{p)(  A(p)  ^  kjp) 

\  +  YCf{p)  \B{p)D{p)  P{p)D{p)  j 


(3.42) 

k  k 

and  /  =  —  - 1  in  the  interval  0  <  /  <  Ymax  -  “  1  •  The  Proof  is  shown  in  Appendix  B. 

k^  k^ 


The  system  Pf^^fip)  in  (3.40)-(3.42)  represents  a  perturbation  of  the  dynamie 
model  (3.24)  due  to  the  presenee  of  the  Low  Pass  Filter  Cp{p).  In  faet,  if  we  let 
Cp{p)  ^1  for  all  p  ,  i.e.,  we  assume  a  filter  with  infinite  bandwidth  the  terms  in  (3.42) 
yield  (  negleet  the  disturbanee  terms  for  eonvenienee ) 


A{p) 


B(p)D(p) 


\+rCp{p) 


and 


kkyo{p)  = 


(1  +  y) 


(\  +  YCp{p))Dip)  Dip) 


This  yields 


y(t)  =  (0^(0 + v{t)) 

Dip) 


(3.43) 


(3.44) 


(3.45) 
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which  is  the  standard  dynamic  operation  in  the  presence  of  parameter  error  6{t)  .  For  the 
input  term,  the  analogy  is  slightly  more  complicated,  but  with  Cp{p)  =  1 ,  simple  algebra 
yields 

=  +  (3.46) 

Now  the  issue  is  how  to  choose  the  filter  C^(5)  so  that  the  transformation 
is  exponentially  stable.  Examination  of  Equation  (3.42)  yields  two  issues: 

1 .  (l  +  /Cp  ( yfy))  *  must  be  stable  for  all  values  of  /  in  the  interval 

0  ^  7  ^  /may  • 

2.  (^7  -  Hyy  ( jco)^  must  be  stable  for  all  allowable  y ,  D{ja)) ,  B{jco) ,  and 
A{jo)) . 

This  leads  to  the  following  two  conditions. 

Condition  1.  Eet  (7»)  be  of  the  form 


CAP)- 


7  ^ 

f  i  ^ 

[p  +  (Oo  j 

C,{p)C,{p) 


(3.47) 


where  L>  N-M  the  relative  degree  of  the  plant.  Also,  let  o,  xug .  Erom  the  Bode  plot, 
shown  in  Eigure  3.3,  it  is  always  possible  to  choose  (U,  and  co^  with  gain  margin  larger 
than  the  parameter  uncertainty  y  . 
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►  CO 


Figure  3.3.  Bode  Plot  of  Filter  Cp{p) . 


For  the  stability  of  (^7  -  (/?)  j  ,  recall  the  expression  for  ( />) 


/^,(rt=(i-c,(p))  1- 


A(p) 


B{p)D{p) 


\  +  YCp{p) 


Since  A{p)  ,  B{p)  ,  and  D{p)  are  all  monic  polynomials  and  B{p)D(p)  is 
Hurwitz  by  assumption,  it  can  be  seen  that 


bm 

p—*OD 


\  Ajp)  '' 

V  B{p)D{p)j 


=  0 


(3.48) 


Since  A{p)  and  [B{p)D{p)')  have  the  same  degree,  this  implies  that 


\  Mp>  ' 

V  B{p)D{p)j 


IS 


low  pass.  At  the  same  time,  with  (5)  being  a  low  pass  fdter  with  unity  DC  gain,  this 
implies  that  (l-C^(/>)]  is  a  high  pass  filter.  Therefore,  the  term 
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1^'  B{p)D{p)j) 


can  be  made  arbitrarily  small  by  inereasing  the 


bandwidth  of  (/>) .  This  is  summarized  by  the  following  eondition; 
Condition  2.  For  every  positive  eonstant  we  ean  always  find  co^  sueh  that 


max< 


JO) 


1-- 


A{jo)) 


B{jG))D{j(o) 


=  <1 


(3.49) 


®  [|7®  +  ®o| 

for  all  CO . 

These  two  eonditions  are  at  the  basis  for  ehoosing  the  filter 

f  V  ^  ^ 

C,{p)=  (3-50) 

\P  +  G)a)y[p  +  Q)^)  j 

given  any  arbitrary  positive  eonstant  <  1 ,  and  ehoosing  cOq  to  satisfy  Condition  2, 
( i.e.,  Equation  (3.49)).  Then  ehoosing  o, ,  so  that  the  gain  margin  is 


Cf  <|l  VMAx\dB  (3.51) 

where  co^  is  the  phase  erossover  frequeney  of  the  filter,  guarantees  that  both  Conditions  1 
and  2  are  satisfied. 

E.  ADAPTIVE  LOOP  FEEDBACK  FUNCTION 

The  parameter  error  term  in  Figure  again  defined  as 

eg{t)  =  Cp{p)(^^ 

This  relates  the  parameter  error  6{t)  with  the  error  eg{t)  and  it  represents  a 
transformation  Hg  ;  ^zi(t)  eg{t)  . 
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Recall  the  key  results  concerning  parameter  estimation  previously  presented  in 
Chapter  III.  C.  In  the  case  of  no  external  disturbances,  the  adaptation  algorithm  yields 
the  following  upper  bounds 

1^1  (0|  11^(0 

and 

d{t)  <//^2(0 

where 

X(t)  =  C,(p)m  (3.54) 

Both  terms  £^{t)  and  are  in  and  their  magnitudes  are  upper  bounded  as 

r,  ^(0)  ' 

j|f,.(r)|jr< -  (3.55) 

The  goal  of  this  is  to  show  that  all  the  bounds  above  (3.52)-(3.55)  yields  an  upper 
of  the  error  e^(t)  due  to  adaptation  as 


\egit)\<  Pit)  sup  \\(piT)\\ 

0<T<t 

(3.56) 

where 

lim  P(t)  =  0 

?— >00 

(3.57) 

and 

r®  V 

\Pit)\<K  — 

(3.58) 

where  is  a  constant  depending  on  the  relative  degree  L,  is  the  bandwidth  of  the 

fdter  Cp  (p) ,  and  is  the  bandwidth  of  the  fdter  (p)  . 

This  is  very  important  since,  going  back  to  Figure  3.2,  equation  (3.56)  shows  that 
the  effect  of  adaptation  e^(t)  is  bounded  by  a  feedback  gain  /]{t) ,  which  not  only  goes  to 
zero  (3.57),  thus  making  the  system  globally  asymptotically  stable,  but  also  its  maximum 


(3.52) 

(3.53) 
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value  (3.58)  depends  on  the  bandwidths  of  the  two  filters  and  it  ean  be  made  arbitrarily 
small.  The  rest  of  this  seetion  is  devoted  to  prove  this  result  and,  not  surprisingly,  is  very 
teehnieal  in  nature. 

The  abounds  (3.56)-(3.58)  will  be  shown  as  sequenee  of  faets. 


Fact  1: 


ee{t)  = 


'csp)'' 


eAt)  = 


kCAp)j 


-V 

< 

V 

0<A<t 


1^1  (0| + p\  1^2 (^)|  Sl  (®£  (^ -  ^))d^  sup  ||^/)(;l)|| 

Q  J 0<A<f 

where  ^^(Oand  determined  by  Equations  (3.33)  and  (3.34),  and 


(3.59) 


L-\ 


Sl  (^) = 


(T-1)! 


^dT  = 


(  L-\  k 

V-t=o  j 


Proof:  First,  rearrange  Equation  (3.39) 


(3.60) 


,{t)  =  CAp)C,{p)(e\t)m\ 


^cap)} 

cap) 


CAp)[cAp)[oAt)m) 


Then  the  error  ean  be  written  as 


eg{t)  = 


C 


ip)[[eAt)CAp)m)+[cAp)[oAm))-[oAt)CAp)m) 


kCAp)j 

If  the  term  in  the  eurly  braeket  is  defined  as  e,(t) ,  this  expression  beeomes 


eg{t)  = 


^C,ip)^ 

KCEiP)j 


(3.61) 


The  term  e,(t)isthe  sum  of  two  terms 


=  +  (0 
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The  first  term  ean  be  written  as 

(0  =  C,{p)[e\t)C,{p)m)  =  C,{p)~e^  (0  (3.62) 

with 

(0  =  e\t)C,{p)m  =  e\t)X{t).  (3.63) 

Additionally,  the  second  term  can  be  expressed  as 

e^(0  =  Co(;7)[Q(;7)(^"(0^(0)-(^"(OQ(;^¥(0)]  =  Co(;?)e^(0  (3.64) 

with 

(0 = Q  {p)[e^  im))  -[o"  it)c,  {p)m\  (3 .65) 

From  equation  (3.33)  we  bound  (3.63)  as 

(0|  =  \e\t)C,{p)m\  =  \o\t)X{t)\  <  I^,(0|||^(0||  (3.66) 

Since  X{t)  =  Cj,{p)(/>{t) ,  the  following  expression  can  be  written: 

dr  sup 

J  0<A<t 

By  the  way  we  choose  the  fdter  Cj,{p)  in  (3.25),  the  impulse  response  (t)  >  0 
for  all  t,  so  the  rightmost  integral  is  C£(0)  which  is  equal  to  one.  Therefore,  X{t)  can  be 
bounded  as  follows: 

|X(t)|  <  sup  ||^z^(/l)||.  (3.67) 

0<yl<f 

Combine  this  with  (3.66)  to  obtain  the  upper  bound 

^|^i(0|||^(0||^K(0|sup||^(2.)||.  (3.68) 

0<x<t 

This  shows  the  bound  of  the  Fact,  in  (3.59). 

For  the  bound  on  (t) ,  expressed  in  Equation  (3.65),  it  is  shown  in  Appendix  C  that 


i  I 

|^(0|  ^  j \ce  {t  -  ^)|  ||^(^)||«?2'  ^  j  \ce  (^) 
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^  c^{T)(l){t -T^dr  dX 


(3.69) 


The  derivative  6’(t)  is  upper  bounded  as  in  (3.53),  so  that  (3.69)  above  beeomes 


eAt)  <^i\\s^{X)\  f  c^{T)dT  dX  sup  ||^(ro) 


(3.70) 


where  again  0^(0  >  0  is  the  impulse  response  of  the  fdter  CAp)  = 


{p  +  g)e) 


- — -  .  Then 


+00  L  L-\ 

j  cAr)dT=  j  = 


(3.71) 


This  shows  the  seeond  inequality  in  (3.59),  and  Faet  1  is  proved. 


Fact  2:  The  error  signal 


QED 


eg{t)=  eXt) 

[cap))' 


(3.72) 


ean  be  upper  bounded  as 


A)  <A0sup||^(r)|| 


(3.73) 


where 


lim  /3{t)  =  0 


/U  \  COq 


\p{t)\<  4  ^  p(0)  ,forallt>0 

I  2)1//  I"  II 


(3.74) 
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A  = 


+^\L-\ 


^tik\ 

Proof:  From  Equation  (3.59) 


2  A 

-lx  r 

e  ax 

j 


•sA  kA  1  [k^+k^y. 

K\kA  ) 


eg{t)  = 


kCAp)j 


eAt) 


with 


e/0  =  e^(0  +  e^(0 

Using  the  bounds  in  Equation  (3.59)  yields 


(0  ^  A(0sup||^(A)|| 


0</l<? 


Where,  using  Schwartz  inequality, 


lAWhj 


=  I  co^e 


-a>„{t-T)  I 


^2  f  +®  2  ^2 


?j(r)<ir<Oo  j  e  J|fj(r)  dr 


This  implies  that 


And,  therefore 


lA(0|s 


(  2 

1  ( 

^(0)  ' 

^-2®o 

0  J 

V 

2//  J 

for  all  t  >  0  .  Eurthermore,  sinee  (t)  e  ,  then 

lim  y0j  (t)  =  0 

?— >00 

Therefore,  Equation  (3.75)  ean  be  rewritten  as 


(3.75) 


(3.76) 


(3.77) 


(3.78) 


(3.79) 
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For  the  second  part  of  eg{t) ,  recall  from  Equation  (3.59) 


with 


^  A(0sup||^(r)|| 


0<r<? 


(3.80) 

(3.81) 


Since  s^it)  e  ,  then 


/->oo 


Again  using  Schwartz  inequality  in  the  convolution  integral  in  (3.81),  we  obtain 
an  upper  bound  for  (t) 


2  f  +-X 


J  V  0 


<A,. 

L  \ 


I  A 

2a>^ 


m 


(3.82) 


Since  Cg(t)  =  Ogg  ,  for  t  >  0 ,  it  is  easy  to  see  that 


A(0  =  Co(0*A(0 


(3.83) 


limy02(O  =  O 


still  yields  the  limit 
and  the  same  upper  bound 

1 

/  +00  '^2 

|A(0|^A  dX  ^\g^{o)^xf  dX 


<A, 


A 


^(0) 


(3.84) 


(3.85) 


Combining  Equations  (3.78)  and  (3.85)  yields  the  result  in  Eact  2. 

QED 

As  noted,  the  dimension  of  the  adaptive  gain,  ju>0  ,  is  inverse  time  or  hme  ' ,  as 
demonstrated  in  Equation  (3.31).  As  a  consequence,  —  in  Equation  (3.74)  is 

Q) 

dimensionless. 

The  bound  on  the  term  e^{t)  can  be  found  in  Equations  (3.73)  and  (3.74)  as 


k/0  ^  AO  sup|Ar)| 

0<r<? 


(3.73) 
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(3.74) 


\fi{t)\<[A,  l^  +  -  ©  ^(0)  ,  foralH>0 

From  Equations  (3.73)  and  (3.74),  e^(t)  can  be  minimized  by  ehoosing  the  adaptive  gain 
//  in  Equation  (3.31)  in  terms  of  the  parameters  cOq  ,  cOj,,  and  degree  L  of  the  fdter 
Cj,  (5)  .  And  then,  by  straightforward  differentiation  the  smallest  upper  bound  is  obtained 

whieh  yields  the  bound  on  the  adaptive  transformation  Hg{p) 

—  sup  ||^/i(r)||  (3.87) 

0)^  ^  0<T<t 

for  some  eonstant  ,  dependent  only  on  the  filter’s  order  and  degree,  L>  N-M  . 

This  result  shows  how  the  upper  bound  on  the  transient  response  e^(t)due  to  adaptation 
depends  in  a  fairly  simple  fashion  on  the  ratio  of  the  bandwidths  and  a>^  of  the  two 

f  CO  ^ 

filters  Cp{s)  and  C^{s)  in  the  eontrol  and  the  estimation  loops.  As  — -  0  the  effeet 

of  adaptation  in  both  input  and  output  signals  tends  to  zero 

In  the  next  seetion,  we  show  this  result  formally,  based  on  the  what  has  been 
shown  in  this  seetion,  summarized  below: 

e„(0  =  ;^e,(0  (3.88) 

Ce(p) 

it  ean  be  shown  that  eg{t)  has  the  same  bound  as  e^{t) 


(3.86) 
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f 


< 


V 


\ee{t)\<  P{t)su^y>{T) 

0<T<t 


m<K 


^  CO.  ^ 


]_ 

4 


\^e) 


lim  P{t)  =  0 

t^oo 


(3.89) 


F.  STABILITY  AND  BOUNDEDNESS  OF  THE  ADAPTIVE  SYSTEM 

In  this  section,  the  concepts  up  to  this  point  are  combined  and  culminate  in  the 
main  result  of  this  dissertation. 

First  define  the  vector  z(t) 


yit) 

u{t) 


(3.90) 


Assuming  zero  disturbances,  from  (3.40)  and  (3.89)  we  obtain 


i  (3.91' 

[  eg{t)^Hg{p)(^^{t)) 

where  is  exponentially  stable,  and  Hg{p)  is  bounded  by  /3{t)  in  Equation 

(3.89).  Also,  recall  that  (/){t)  is  a  vector  of  filtered  input  and  output  signals,  which  can  be 
written  as 

m  =  ^{p)^{t)  (3.92) 

where  0(/?)  is  exponentially  stable  transformation  in  (3.17)  and  (3.18). 

The  whole  system,  obtained  by  combining  Equations  (3.91)  and  (3.92),  is  shown 
in  Eigure  3.4. 
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Figure  3.4.  Dynamics  Model  of  the  Adaptive  System  with  No  Disturbance. 

Since,  ^  0  (Equation  (3.89)),  and  and  ^(p)  are  exponentially 

stable,  the  closed  loop  system  is  exponentially  stable  and  all  signals  are  bounded. 

Also,  if  the  reference  input  and  output  signal  are  defined  as 

z,,,{t)  =  H^^^{p)CAp)v{t)  (3.93) 
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From  these,  it  is  easy  to  see  that 


ML  ^  |®(;^)L  (IML  +||v|L)  (3.99) 


whieh  yields 


(3.100) 


K=H,.,{p)\H,{p)Ui>{pl 


(3.101) 


Inequality  (3. 100)  is  valid  only  when  the  gain  K  is  smaller  than  one  as  0  <  <  1 . 

From  the  results  above,  reeall 


\Hg{p) ^  =  /3{t)  <k^ 


(3.102) 


where  depends  on  the  filter  order  only.  Therefore,  it  is  observed  that  as  — ^  q  ,  the 


overall  loop  gain  K  goes  to  zero  and,  therefore,  the  error  terms 


z(0  = 


yit)- Prefix) 
u{t)-uAt) 


(3.103) 


are  such  that 


lim  sup  y(0-T../(0  =0 


_Q — >00  0<?<+oo 
COe 


(3.104) 


lim  sup  M(t)-M^^^(t)  =0. 

^ _ >00  0<f<+oo 


(3.105) 


In  other  words,  the  maximum  value  of  the  transient  response  due  to  the  adaptive 
controller  tends  to  become  arbitrary  small. 
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In  the  next  seetion,  we  show  that  eomputer  simulations  support  the  results 
presented  above. 

G.  EXAMPLES  OF  APPLICATION 

1,  LI  Adaptive  Control  in  Ideal  Case  (No  Disturbance) 

Consider  an  unstable  plant 

yif)  =  K-^u{t)  (3.106) 

p-1 

where  1  <k,<  4 ,  and  the  referenee  model 


p  +  2 


(3.107) 


Then 


B(p)  =  l  ,A(p)  =  p-l,  D(p)  =  p  +  2. 


(3.108) 


First,  the  true  eoeffieients  are  determined,  assuming  =3  as  the  true  value  of  the 
parameter.  From  A(p)  and  B(p) ,  the  related  spaee  state  equation  ean  be  determined  and 
is  expressed  as  follow: 


x{t)  =  A^x{t)  +  BpU{t) 
y{t)  =  CpX{t) 


where 


Following  the  steps  mentioned  above,  the  elosed  loop  system  is  obtained 

x{t)  =  {Ap  -BpLp)x{t)  +  BpV{t) 

y{t)  =  CpX{t) 


where  eontroller  gain  =3  is  derived  from  Equation  (3.9). 
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Let  the  observer  be 


x{t)  =  ApX{t)  +  BpU{t)  +  Kp  {y{t)  -  y{t)) 

y{t)  =  CpX{t) 


(3.111) 


where  P{s)  =  p +  2  and  the  observer  gain  =  1  is  derived  from  Equation  (3.10). 

From  Equations  (3.15),  (3.16),  (3.22),  and  (3.23),  the  system  ean  be  represented 


as 


(3.112) 


where 


=[2  3  4.5f 


(l){t)  = 


u{t) 


s  +  2 


u{t) 


s  +  2 


yit) 


(3.113) 


Now,  with  kQ=3,  a  suffieiently  large  value  for  co^  must  be  ehosen.  Additionally, 
aeeording  to  Condition  1  and  Equation  (3.49),  the  magnitude  of  the  response  of 


\  Ap)  ^ 

V  B{p)D{p) 


( 


1-- 


A{p) 


must  be  less  than  1 .  The  magnitude  of  the  Bode  Plot  of 


is  shown  in  Figure  3.5  and  it  is  smaller  than  one  (negative  dB’s)  for  all 


B{p)D{p) 

frequeneies  larger  than  than  2.5  rad/see.  This  ean  be  used  as  a  guideline  to  determine  the 
bandwidth  co^  of  the  fdter  C^(5) . 
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bode  of  (1-(A(s)/(B(s)D(s))) 


Figure  3.5. 


The  Bode  Plot  of  1 

V 


A{p) 

B{p)D{p)^ 


With  a  choice  of  Oq  =  5  ,  the  LPF  becomes  Cp  {p)  = 
must  verify  that  Condition  2  is  satisfied.  Recall  this  condition 


5  +  5 


and  we 


where 


max  max  \h  ( jco)\  <  1 

®  o<r<rMAx '  ' 


h„{p)={i-cAp)) 


A{p) 


B{p)D{p) 


^+yC,[p) 


y 


For  a  number  of  different  realizations  of  y,  is  plotted,  as  shown  in 

Figure  3.6.  The  simulation  results  conform  to  the  requirement  of  Condition  2. 
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bode  of  Hyy{s) 


Figure  3.6.  The  Bode  Plot  of  H^{p)  for  100  Plant  Realizations. 


Next,  if  cOj,  is  chosen  such  that  cOj,>cOq,  and  the  degree  L  of  C,,{p)  is  guaranteed  to 
meet  L  >  N-M ,  then  C^{p)D{p)  will  be  proper.  In  this  example,  let  cOj,  =100  rad/sec, 
L  =  \  =  M -N ,  and,  therefore. 


cAp)  = 


a>^ 


100 


ip  +  (i>E)  (;?  +  100) 


(3.114) 


Since  this  case  is  implemented  with  no  disturbances,  W£(t)  =  0  .  Figure  3.7  and  Figure 
3.8  depict  simulation  results  of  the  case  without  disturbances. 
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1 


input  u(t) 


0.8  ^ 

0.6  -  - 

0.4  -  - 

0.2f  I 

-0.2  -  '  '  '  -I 

-0.4  -  I 

-0.6  - 

-0.8  -  I 

_-|  _ I _ I _ ^ _ I _ ^ _ I _ I _ ^ ^ _ 

0  5  10  15  20  25  30  35  40  45  50 

sec 

Figure  3.7.  The  Simulation  Result  of  the  Control  Input. 


reference  output  and  plant  output 


Figure  3.8.  The  Simulation  Result  of  the  Referenee  Output  and  Plant  Output. 
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2,  LI  Adaptive  Control  with  Bounded  Output  Disturbance 

Consider  a  system 

p-1 

where  0.5  <k,<  3  ,  and  the  referenee  model  is 

p  +  \ 

where 

D(p)  =  p  +  \. 

Then 

B{p)  =  \  ,A{p)  =  p-\  ,D{p)^p  +  \  (3.115) 

Assuming  =  2  as  the  true  value,  the  eontroller  ean  be  expressed  as 

u{t) 

0.5u(t)  =  -^kQ-k^  A:^]  +v{t)  =  -^  ^{t)  +  v{t)  (3.116) 

[l_ 


where 


’^^=[0.5  2  6f 

(i){ty  =  u{t)  m(o  t(o 
p +  5  p +  5 


(3.117) 


In  aeeordanee  with  Condition  1  and  Condition  2,  terms  are  ehosen  as  (0^=5 , 
(y,  =  10  ,  =  100 ,  and  L  =  2>N-M  =  1.  As  in  the  previous  example,  Figure  3.9  and 

Figure  3.10  show  that  the  stability  eonditions  are  satisfied. 

With  appropriate  substitutions  (/>)  and  Cj,{p)  ean  be  expressed  as: 
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CM  = 


COr, 


yp  +  Q)J 


CO, 


(p  +  co^) 


f  5  ) 

f  ^ 

10" 

Ce(p)  = 


co^ 


(P  +  COe) 


100 

(;?  +  100)' 


bode  of  (1-(A(s)/(B(s)D(s))) 


Frequency  (rad/sec) 


Figure  3.9. 


f 

The  Bode  Plot  of  1 

V 


A{p) 

B(p)D{p)^ 


bode  of  Hyy{s) 


Figure  3.10.  The  Bode  Plot  of  Hyy{p)  for  100  Plant  Realizations. 
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Figure  3.11  and  Figure  3.12  show  the  results  of  this  simulation. 


input  u(t) 


Figure  3.11.  The  Simulation  Result  of  the  Control  Input. 


Figure  3.12  The  Simulation  Result  of  the  Reference  Output  and  Plant  Output. 
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3,  Example  with  Sensor  Noise 

Consider  a  plant  with  transfer  funetion 


A{p)  p  +  a 

where  the  parameters  are  within  known  bounds 

-2  <  a  <  0 ,  0.5  <  kn  <  3 


Let  the  referenee  model  be 


1 


1 


Dip)  p  +  l 

First,  the  LPF,  C^ip)  is  chosen  and  is  of  the  form 


cAp)  = 


f  \ 

(  i  ^ 

0)^ 

yp  +  O)^  j 

{{p  +  (D,f  ] 

=  C,ip)Qip). 


Again,  we  choose  L  =  2>N-M=\  .  From  the  bound  on 


0  <  7  <  7^^  =  —  - 1  =  — — 1  =  5  and  the  Bode  plot  of 

0.5 


1-- 


Aip) 


Bip)Dip) 


Figure  3.13  for  a  number  of  random  realizations  of  the  plant. 


(3.118) 

(3.119) 

(3.120) 

(3.121) 

,  we  obtain 
is  shown  in 
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Figure  3.14.  Frequency  Response  of  H^y{p)  for  100  Plant  Realizations. 
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It  can  easily  be  seen  that  for  ty  >  3  rad/see  all  realizations  have  magnitude  smaller  than 
one  that  required  by  Condition  1.  In  this  example,  (yg  =5  rad/see  is  ehosen  so  that 

Condition  1  holds  for  all  realizations,  and  (y,  =10  rad/see  so  that  the  gain  margin  of 
Cp{p)  is  larger  than  lOlogjg  =  \AdB .  This  idea  ean  be  applied  to  the  eontroller 
design. 

From  Figure  3.14,  the  magnitudes  of  all  realizations  are  found  to  be 

smaller  than  one  required  to  satisfy  Condition  2.  Both  of  these  results  guarantee  the 
stability  of  in  Equation  (3.40)  for  all  realizations  of  the  plant. 


The  responses  of  the  inputs  and  output  signals  for  this  adaptive  eontrol  system  for 
two  different  values  of  co^-  are  shown  in  Figure  3.15  and  Figure  3.16.  To  ensure  a  more 
realistie  experiment,  bandlimited  random  observation  noise  is  added.  The  simulation 
results  are  shown  in  Figures  3.17  and  Figure  3.18.  Notiee  that,  as  predieted  by  the 
theoretieal  arguments  in  the  previous  seetion,  the  peak  of  the  transient  response  deereases 
as  the  bandwidth  cOj,  inereases.  The  relationship  between  co^,  cOj,,  and  the  transient 


performanee  is  evident  in  Equation  (3.87). 


''  CO,  A 


V  J 


sup  (/){t) 


a<T<t 


(3.122) 


In  the  next  seetion,  the  upper  bounds  of  parameters  and  transfer  funetions  is  diseussed. 
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Figure  3.15.  Frequency  Response  to  Square  Wave  for  cOj,  =10  rad/sec. 
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Figure  3.16  Frequency  Response  to  Square  Wave  for  =100  rad/sec. 
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input  with  wE  =  10  rad/sec  output  with  wE  =  10  rad/sec 


Figure  3.17.  Frequency  Response  to  Square  Wave  for  =  10  rad/sec  with 
Bandlimited  Random  Observation  Noise. 


input  with  wE  =  100  rad/sec  output  with  wE  =  100  rad/sec 


Figure  3.18.  Frequency  Response  to  Square  Wave  for  co^  =  100  rad/sec  with 
Bandlimited  Random  Observation  Noise. 
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4,  Example  with  Modeling  Error 

In  this  case,  the  adaptive  controller  presented  in  the  previous  section  is  used,  but 
the  plant  has  unmodeled  dynamics  as  represented  by 


yit)  =  ko 


1  200 
p-a  /7^+30/?  +  200 


u{t) 


(3.123) 


In  other  words,  the  system  is  treated  as  a  first  order  system,  while  in  reality  it  has  two 
extra  poles  affecting  the  high  frequency  only.  It  is  well  known  that  in  this  situation  the 
adaptive  controller  can  easily  become  unstable  and  sensitive  to  external  disturbances. 
However,  the  proposed  controller  seems  to  retain  stability  and  traeking  properties. 
Simulation  results  are  shown  in  Figure  3.19  and  Figure  3.20. 
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Figure  3.20. 
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IV.  LI  ADAPTIVE  CONTROL  FOR  NON-MINIMUM  PHASE 

SYSTEMS 

A,  MOTIVATION  AND  BACKGROUND 

The  direct  Adaptive  Control  technique  presented  in  the  previous  chapters 
assumes  that  the  system  to  be  controlled  is  minimum  phase,  in  the  sense  that  the  zeros  of 
its  transfer  function  have  to  be  on  the  left  half  side  of  the  s-plane. 

It  is  well  known  that  non-minimum  phase  zeros  present  very  important 
limitations  on  the  performance  of  control  systems.  Most  important  is  the  fact  that,  while 
poles  can  be  moved  by  state  feedback,  zeros  can  be  replaced  only  by  pole-zero 
cancellation.  Since  unstable  pole-zero  cancellation  leads  to  unstable,  uncontrollable,  or 
unobservable  modes,  this  operation  is  clearly  not  feasible. 

Non-minimum  phase  systems  are  very  important  in  some  applications  and 
therefore  cannot  be  ignored.  In  particular,  systems  with  flexible  appendages  most  likely 
have  zeros  on  the  imaginary  axis  and  are  therefore  non-minimum  phase  [89]. 

In  this  chapter,  the  adaptive  control  of  this  class  of  systems  is  discussed.  The 
non-minimum  phase  limitation  is  avoided  by  controlling  an  augmented  plant  obtained  by 
the  addition  of  a  system  in  parallel  to  the  plant  to  be  controlled.  This  augmenting 
network  has  to  be  designed  such  that  its  output  is  zero  in  steady  state,  and  the  overall 
combination  is  minimum  phase. 

Although  the  control  presented  for  minimum  phase  systems  is  restrictive,  it  does 
exhibit  desired  tracking  in  the  presence  of  wide  plant  uncertainties.  Several  computer 
simulations  show  the  effectiveness  of  such  a  controller.  Furthermore,  implementation  on 
an  experimental  flexible  mechanical  arm,  with  uncertain  dynamics,  shows  that  this 
approach  can  be  a  viable  technique  for  the  control  of  this  class  of  systems. 


B,  FUNDAMENTAL  PRINCIPLE  AND  ANALYSIS  OF  THE  MODIFIED  LI 
ADAPTIVE  CONTROL 

Along  similar  lines  as  in  the  previous  chapters,  consider  a  Single  Input  Single 
Output  (SISO)  system  with  dynamics 
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y{t)  =  k,^u{t). 
A(p) 


(4.1) 


with  possible  non  minimum  phase  zeros.  When  this  is  the  case,  as  shown  in  Figure  4.1, 
we  augment  it  by  a  parallel  system  so  that  the  combination  system  becomes 


m= 


V 


eAp)'" 

Dl{p)y 


u{t)  = 


Ap)E,  ip)  +  KB{p)D^  ip) 
A{p)D,{p) 


(4.2) 


Figure  4.1. 


New  Closed  Loop  Dynamics  with  a  Designed  PID 


E,{p) 

E>i{p) 


Now,  the  problem  is  to  choose  the  augmenting  transfer  function 


Eiip) 

E>i{p) 


to  satisfy  the 


following  requirements: 

a. )  the  numerator  of  the  combined  system  {^A{p)E^{p)  +  k^^B{p)D^{p)^  is  Hurwitz; 

b. )  in  steady  state  the  output  of  the  augmenting  system  goes  to  zero. 

To  satisfy  Requirement  b.)  when  the  input  is  a  constant,  Ei{p)  is  expressed  as 


E,{p)  ^  P'El  {p) 


(4.3) 
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so  that  the  derivative  aetion  (the  term  “  p  “)  makes  the  steady  state  value  zero  when  the 
input  is  a  constant.  To  satisfy  Requirement  a.),  notice  that  the  numerator  of  Equation 
(4.2)  is  the  characteristic  polynomial  of  the  system  shown  in  Figure  4.2. 


u{t) 


E>,{P) 

- ^ 

KB{p) 

i 

L 

E,{p) 

A{p) 

- P* 

t(0 


Figure  4.2.  The  Other  Expression  of  the  New  Closed  Foop  Dynamics  with  a  Designed 

PID. 


Therefore,  the  augmenting  network 


Eiip) 

D,{p) 


must  be  such  that  its  inverse 


Dxip) 

E,{p) 


yields  a 


stable  closed  loop  system. 


Furthermore,  since  it  is  desired  that 


eAp) 

E>,{p) 


is  strictly  proper. 


E>,{p) 

Exip) 


will  not  be  proper. 


A  possible  choice  is  that 


E>i{p) 

E,{p) 


is  a  PID  controller,  with  the  integral  action  to 


satisfy  Equation  (4.3),  and  the  derivative  action  to  make  its  inverse  proper.  For  example, 
let  a  nominal  plant  be  expressed  as 


B{p)  _  /+!  {P  +  V){P-V) 

A(p)  °  p'^ +2p^  ^  p^i^p  +  \A\A2j^i^p-\A\A2j^ 


Observe  that  there  are  two  zeros  located  on  the  imaginary  axis  and  therefore,  it  is  non¬ 
minimum  phase. 
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In  order  to  make  the  system  minimum  phase,  a  PID  compensator  must  be 

E,{p) 

designed  so  that  the  system  in  Figure  4.2  is  stable.  The  augmented  minimum  phase 
system  is  shown  in  Equation  (4.2). 

From  Figure  4.1,  the  new  output  is  shown  to  be 


yit)  =  yit)  y^ui.)  = 


D,{p) 


A{p)D,{p) 


(4.5) 


Based  on  the  requirements  listed  above,  the  degree  of  E^{p)  has  to  be  less  than£)j(/»). 
Further,  according  to  the  final  value  theorem  [89-90],  the  numerator  E^{p)  must  have 
roots  at  its  “zeros“.  Thus  define 


Ey{p)  =  P'E^ip) 


Z)j {p)  =  UqP^  +  a^p‘^  '  H - ha 


and  choose 


Ei{p)  =  p 


Dy{p)  =  p^  +  p  +  Q.\ 


Then  the  transfer  function  of  the  augmented  system  becomes 

BB(p)  _  A{p)E^ {p)  +  B{p)D^ {p)  _  p^  +  p"^  +  3p^  +  \.\p^  +  p  +  i).\ 


AA(p)  A(p)mp) 


p^+p'+2.1p^+2p^+0.2p^ 


with  zeros  at 


(4.6) 

(4.7) 

(4.8) 

(4.9) 


(4.10) 


p  =  [-0.3347  +  1.5021i  -0.3347  -  1.50211  -0.1107  +  0.61141  -0.1107  -  0.61141  -0.1094] 


(4.11) 


Which  are  strictly  located  in  the  FHS  s-plane.  This  leads  to  the  augmented  system 
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(4.12) 


m= 


BB{p) 

AA{p) 


u{t)  = 


p^  +  p*  +3p'  +\.\p^  +  p  +  OA 

+2.1/ +2/ +0.2/  ^ 


Now,  the  reference  model  is  chosen  as  in  Equation  (3.4): 

p  +  \ 


(4.13) 


where  D(p)  =  p  +  1  is  a  Hurwitz  polynomial  representing  the  desired  closed  loop 
dynamics.  The  degree  of  D(p)  is  N-M ,  which  is  equal  to  the  relative  degree  of  the 
augmented  plant,  and  v(t)  is  an  arbitrary  external  bounded  input. 


Designing  the  standard  adaptive  controller  for  the  augmented  minimum  phase 
system  follows.  In  particular,  let  the  filter  be  chosen  as 


cAp)  = 


f  \ 

(  2  ^ 

r  10  ^ 

f  0  ^ 

10" 

[p  +  coj 

[;7  +  loJ 

j 

(4.14) 


It  is  possible  to  verify  that  this  guarantees  stability  of  the  reference  model. 


max  max  H  ( jco)  <  1 

(9  o<r<rMAx 


(4.15) 


where  H^[p)  is  as  defined  in  Chapter  III.  Simulation  results  shown  in  Figures  4. 3-4. 5 
show  the  stability  and  asymptotic  tracking  of  the  desired  trajectory. 
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modified  plant  output  and  original  plant  output 


Figure  4.3.  Modified  Plant  Output  and  Original  Plant  Output. 
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Figure  4.4.  Simulation  Result  of  Modified  LI  Adaptive  Control  System. 
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Figure  4.5.  Modified  Plant  Output  Tracks  the  Reference  Output. 


C.  SIMULATIONS  OF  MODIFIED  LI  ADAPTIVE  CONTROL  WITH  NON¬ 
MINIMUM  PHASE  SYSTEM 

In  order  to  verify  the  robustness  and  stability  of  the  modified  LI  adaptive 
controller  with  PID  element,  a  modified  controller  is  adapted  to  perform  under  realistic 
conditions. 


1,  Modified  LI  Adaptive  Control  System  with  Bounded  Output 
Disturbance 

Consider  a  system 


y{t)  =  K 


Bjp) 

A(p) 


u{t)  = 


0.06023/ +0.13378 
/ +2.57184/ 


u{t) 


(4.16) 
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where  0.5  <k,<3  ,  assuming  A:q  =  1  and  the  referenee  output 


p  +  l 


Then 

5(;7)  =  0.06023/ +0.13378,  A{p)  =  +  2.5mA p\  D{s)  =  p  +  \  (4.17) 

This  plant  has  two  zeros,  one  loeated  at  (0  +  1 .49041)  and  the  other  at  (0  -  1 .49041),  both 
loeated  on  the  imaginary  axis,  and  therefore  the  system  is  non-minimum  phase. 

E  (  £?) 

Next,  a  proper  PID  eontrol  — — -  is  defined  with 

mp) 


and 


E^(p)  =  p 


(4.18) 


Dy{p)  =  \Qp^  +5p  +  Q.\ 


(4.19) 


It  is  easy  to  verify  that  this  eontroller 


E,{p) 

E>i{p) 


stabilizes  the  system. 


After  the  system  is 


augmented  with 


E,{p) 

E>i{p) 


the  response  is 


,  BB{p)  / +0.6023/ +2.8730/ +1.3438/ +0.6689;7  +  0.0134 
y{t)  =  ^0 - —u{t)  =  E - - - 1— - - i— - i— - i— - u(t). 

AA(p)  /+0.5/+2.5818/  +  1.2859/+0.0257/ 

(4.20) 

Now,  the  referenee  model  is  ehosen  as 

=  (4.21) 

P  +  5 

where  D{p)  =  (^p  +  5)  is  a  Hurwitz  polynomial  representing  the  desired  elosed  loop 
dynamies.  The  degree  of  D{p)  is  N-M ,  whieh  equals  the  relative  degree  of  the  new 
plant,  and  v(t)  is  an  arbitrary  external  bounded  input. 


68 


Furthermore,  in  order  to  satisfy  the  stability  eonditions  of  Chapter  III,  Oq  =  30 


rad/see,  <»,  =  30  rad/see,  and  =  50  rad/see  are  ehosen.  Then  the  LPF  beeomes 


C,{p)  = 


f  \ 

(  2  ^ 

nij 

r  30  ^ 

f  9  A 

30" 

yp  +  a)J 

l;?+3oJ 

IL+sofJ 

(4.22) 


and 


cAp)  = 


{p  +  COEf 


50" 

(;7  +  50)"' 


We  ean  verify  that,  with  this  LPF, 


max  max  (  ico)  <  1 

m  0<y<ruAx 


(4.23) 


(4.24) 


and  the  transfer  funetions  '.H^i^p^,  ean  be  obtained 

from  Equation  (3.42)  as  before. 

Figures  4. 6-4. 9  show  the  results  of  this  simulation.  From  Figure  4.6,  it  ean  be 
seen  that  the  augmented  plant  is  stabilized.  In  the  presenee  of  a  bounded  output 
disturbanee,  the  settling  time  is  a  little  longer  eompared  to  the  example  with  no 
disturbanee  in  last  seetion.  However,  due  to  the  effeet  of  the  additional  PID  element,  the 
original  output  tracks  the  modified  plant  output  and  it  is  still  stable,  despite  the  output 
disturbance.  In  addition,  from  Figures  4.7  and  4.8,  it  is  obvious  that  this  modified  LI 
adaptive  control  system  satisfies  the  specified  requirements  of  this  experiment.  It  is 
evident  that  there  is  no  large  control  action  or  poor  transient  performance  during  the 
duration  of  execution,  and  the  modified  plant  output  did  track  the  reference  output. 
Figure  4.9,  shows  that  the  adaptive  error  is  small  and  bounded,  which  satisfies  the 
requirement  of  a  stable  adaptive  control  system. 
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Figure  4.6.  Modified  Plant  Output  and  Original  Plant  Output  with  Bounded  Output 

Disturbance. 
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Figure  4.7.  Simulation  Result  of  Modified  LI  Adaptive  Control  System  with  Bounded 

Output  Disturbance. 
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Figure  4.8.  Modified  Plant  Output  Tracks  the  Reference  Output  with  Bounded  Output 

Disturbance. 


reference  output  and  modified  plant  output 

T  I  ^  I  I 

- reference  output 

- modified  plant  output 


0.1 

0.05 

0 

-0.05 

-0.1 


input  u(t)  error 


300 

sec 

tracking  error 


600 


100 


200 


300 


400 


500 


600 


Figure  4.9.  The  Adaptive  Error  of  Modified  Plant  Output  with  Bounded  Output 

Disturbance. 
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2,  Modified  LI  Adaptive  Control  System  in  the  Presence  of  Unmodeled 
Dynamics 

In  this  case,  the  same  adaptive  controller  is  used  as  in  the  previous  example,  but 
the  plant  has  unmodeled  dynamies  as  represented  by 


y{t)  =  K 


0.06023/ +0.13378  300 

/ +2.57184/  /+10;7  +  30o“^ 


(4.25) 


where  =  0. 1  ,  /  =  0.9  *  =  0.09  <  <  1 . 1  *  =  0. 1 1  is  assumed.  In  addition, 

from  Equation  (4.25),  it  has  two  extra  poles  affecting  the  high  frequencies  only.  In 
general,  in  this  situation,  the  adaptive  controller  can  easily  beeome  unstable  and  sensitive 
to  external  disturbanees.  Also,  the  plant  has  two  zeros,  one  loeated  at  (0  +  1 .4904i)  and 
the  other  at  (0  -  1.4904i),  both  on  the  imaginary  axis  and  therefore  the  system  is  non¬ 
minimum  phase. 


Next,  the  proper  control  is  defined  as 


Eiip) 


,  the  inverse  of  a  PID  eontrol,  where 


Eiip)  =  p 


(4.26) 


and 


D^(p)  =  10j!7^+5/>  +  0.1. 


(4.27) 


E  ( 

It  is  easy  to  verify  that  this  controller  — — -  stabilizes  the  system.  After  augmenting  the 


E>i(p) 


.  ,  E.(p)  , 

system  with  ,  the  response  is 

A(/ 


,  ,  BBip)  ,  / +0.6023/ +  2.8730/ +1.3438/ +0.6689;7  +  0.0134 

=  K - 6  .  ^  .  s  .  .0.0  4  .  .  _ _ _  7 - 


AA(p) 


/  +0.5/  +2.5818/  +1.2859/  +0.0257/ 


(4.28) 


Now,  the  referenee  model  is  ehosen  as 
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(4.29) 


p  +  5 

with 

Dip)  =  p  +  5  (4.30) 

where  Dip)  =  (^p  +  5)  is  a  Hurwitz  polynomial  representing  the  desired  elosed  loop 
dynamies.  The  degree  of  Dip)  is  N-M ,  whieh  equals  the  relative  degree  of  the  new 
plant,  and  v(t)  is  an  arbitrary  external  bounded  input. 

Furthermore,  in  order  to  satisfy  the  stability  eonditions  of  Chapter  III,  =  30 
rad/see,  o,  =  30  rad/see,  and  =  50  rad/see  are  ehosen.  Then  the  LPF  beeomes 


cAp)  = 


r  \ 

f  2  ^ 

^  30  ^ 

f  0  ^ 

30" 

[p  +  «}o  j 

{{p  +  0),f  J 

^;?  +  3oJ 

iL+sofj 

(4.31) 


and 


cAp)  = 


ip  +  (oA 


50" 

Xp^ 


(4.32) 


It  is  easy  to  verify  that,  with  this  LPF, 


max  max  H  ijco)  <\  (4.33) 

w  0<y<r„^x 


and  the  transfer  funetions  '.H^i^p)  ,  H ^^i^p)  ,  and  can  be  easily 

obtained  from  Equation  (3.42). 

Figures  4.10-4.12  show  the  simulation  results  of  this  experiment.  From  Figure 
4.10  and  Figure  4.11  it  is  obvious  that  this  modified  LI  adaptive  eontrol  system  satisfies 
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the  specified  requirement.  It  is  evident  that  there  is  no  large  control  action  or  poor 
transient  performance  during  the  duration  of  execution,  and  the  modified  plant  output  did 
track  the  reference  output.  In  Figure  4.12,  the  adaptive  error  is  shown  and  is  small  and 
bounded,  which  satisfies  the  requirement  of  a  stable  adaptive  control  system. 
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Figure  4.10.  Simulation  Result  of  Modified  LI  Adaptive  Control  System  with 
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Figure  4.11.  Modified  Plant  Output  Tracks  the  Reference  Output  with  Unmodeled 
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D.  APPLICATION  TO  A  SATELLITEWITH  A  FLEXIBLE  ARM 


In  this  section,  the  problem  of  applying  the  proposed  adaptive  controller  to  the 
control  of  a  flexible  arm  is  addressed.  This  is  part  of  the  equipment  available  in  the 
Satellite  Research  and  Design  Center  (SRDC)  at  the  Naval  Postgraduate  School.  The 
SRDC  aims  to  recreate  realistic  space  conditions  in  a  laboratory  environment. 

The  particular  equipment  used  in  this  experiment  is  the  Flexible  Spacecraft 
Simulator  (FSS)  shown  in  Figure  4.13.  The  FSS  simulates  a  flexible  arm  moving  in  the 
plane. 


Figure  4.13.  Flexible  Spacecraft  Simulator  (FSS). 
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It  was  developed  to  experimentally  verify  the  effeetiveness,  robustness  and  stability  of 
designed  eontrol  laws.  It  is  eomprised  of  a  rigid  eentral  body  and  a  flexible  appendage. 
The  rigid  eentral  body  represents  the  main  body  of  the  spaeeeraft,  and  the  flexible 
appendage  represents  a  flexible  antenna  support  strueture  [91].  In  addition,  the  flexible 
appendage  is  eomposed  of  a  based  beam  eantilevered  to  the  main  body  and  tip  beam 
rigidly  eonneeted  to  the  base  beam  at  a  right  angle.  Five  air  pads  on  a  granite  table 
support  the  FSS  and  its  flexible  appendage  to  minimize  frietion  during  motion. 

The  flexible  model  is  similar  to  the  mass-spring-damping  model.  Based  on 
experimental,  real  data  from  FSS,  the  equations  of  motion  of  the  flexible  model  ean  be 
developed  and  are  written  as 


M  {x}  -t-if  |x}  =  B^u 


(4.34) 


and 


M  = 


19.2253 

-1.6190 


-1.6190 

1 


,K  = 

“0 

0 

.  A  = 

“O' 

0 

2.2212 

1 

,  and  X  = 


0 


where  0  is  the  angle  of  the  spaeeeraft  hub  in  radians  and  g  is  a  state  variable 
representing  appendage  deformation  using  modal  eoordinates.  From  Equation  (4.34),  the 
state  spaee  dynamies  are  derived  and  are  written  as 


where 


[  x(t)  =  ApX(t)  +  BpUit) 

1  y(t)  =  CpX(t) 


0  0  10 
0  0  0  1 
0  -2.166  0  0 
0  -2.5718  0  0 
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0 

0 

0.06023 

0.09751 


and 


C,=[l  0  0  0] 

where  y{t)  =  6{t)  is  the  measured  angular  position. 


With  appropriate  substitutions,  the  transfer  function  of  the  plant  becomes 


B{p)  0.06023/ +0.13378 
/O  =  _ _  ,  u{t) 


A{p) 


/ +  2.57184/ 


(4.35) 


It  can  be  seen  that,  this  flexible  model  has  two  non-minimum  phase  zeros  located  on  the 
imaginary  axis.  The  root  locus  plot  of  this  plant  is  shown  in  Figure  4.14. 
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Figure  4.14.  The  Root  Locus  Plot  of  the  Non-minimum  System 
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It  can  also  be  verified  that  the  PID  eontroller  with  transfer  funetion 


D^{p)  _10/>^+5j!7  +  0.1 

E,{p)  p 


(4.36) 


stabilizes  the  system.  Also,  it  provides  suffieient  margin  so  that  parameter  perturbations 
do  not  affect  the  stability.  Then,  the  nominal  augmented  system  is 


m= 


B{p)  ,  E^{p) 


-  + 


Ap)  E)^{p), 


u{t) 


(4.37) 


and  beeomes 


_  ^  BB(p)  _pA  0.6023 /  +  2.8730y  + 1 .3438;?^  +  0.6689 p  +  0.0134 

A  \  /  0  A  A  /  \  ^  .  r\  r  5)  .  r\  r  c\  ■%  c\  A  .■%  r\  c\  r  r\  7)  .  r\  r\r\  r  ^  2  ^  ' 


AA(p) 


/  +0.5/  +2.5818/  +1.2859/  +0.0257/ 


(4.38) 


where  kg  =1  and  A:o  =  0.1 . 


Sinee  the  relative  degree  is  A-M  =  1 ,  the  referenee  model  is  ehosen  as 

y,(0  =  /v(f) 
p  +  5 


(4.39) 


with  v(t)  representing  an  arbitrary  external  input.  It  ean  be  easily  verified  that  the 
augmented  system  is  minimum  phase,  and  ean  be  adaptively  eontrolled  to  follow  the 
trajeetory  of  the  referenee  model  of  Equation  (4.39). 

Proeeeding  as  in  the  previous  ehapter,  we  ehoose  the  fdters  Cj^(s)  and  C^(5)  for 
the  purposes  of  eontrol  and  adaptation.  With  cOg  =  30  rad/see,  co^  =  30  rad/see,  and 
CO =50  rad/see  ehosen,  the  EPF  beeomes 


cAp)  = 


c  \ 

(Og 

(  2  ^ 

r  30  ^ 

f  9  A 

30' 

yP  +  COg) 

l;?+3oJ 

l(/>+3o)y 

(4.40) 


and 
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L 


(4.41) 


cAp)  = 


co^ 


{P  +  COe) 


50^ 

(;?  +  50)^' 


From  the  plot  of  |//_^(70)|  in  Figure  4.15,  it  can  be  see  that  the  magnitude  is 
always  smaller  than  “one”,  thus  guaranteeing  stability  of  the  transfer  functions  :  , 

Hyo  {p) ,  {p) ,  and  {p)  in  Equation  (3.42). 


abs(Hyy(jw )) 


Figure  4.15.  The  Plot  of  the  |//^(7(»)| 


Figures  4.16  and  4.17  show  the  results  of  this  experiment.  Figure  4.16  shows  the 
trajectories  of  the  reference  model,  original  model,  and  the  modified  plant.  Figure  4.17 
shows  the  results  of  control  input,  reference  output,  modified  plant  output,  and  original 
plant  output,  respectively.  From  these  figures,  it  can  be  seen  that  the  modified  plant  with 
PID  control  results  in  the  plant  becoming  stable.  Comparing  the  figure  of  the  reference 
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output,  it  is  evident  that  the  modified  output  after  adding  the  PID  does  track  the  trajectory 
of  the  reference  model.  In  addition,  no  large  control  action  or  poor  transient  performance 
is  exhibited  during  the  entire  execution  period.  This  result  validates  that  the  modified 
adaptive  controller  with  PID  design  does  solve  the  flexible  problems  in  non-minimum 
phase  systems.  This  is  the  reason  that  such  a  controller  is  implemented  in  this 
experimental  adaptive  system. 


Figure  4.16.  Trajectories  of  the  Reference  Model,  Original  Plant,  and  Modified  Plant 
with  Modified  LI  Adaptive  Controller  and  FSS. 
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V.  CONCLUSIONS 


A,  SUMMARY  AND  CONTRIBUTIONS 

In  this  dissertation,  an  Adaptive  Control  system  has  been  redesigned  for  Single 
Input  Single  Output  (SISO)  systems  to  improve  both  transient  response  and  robustness  to 
external  disturbanees  and  unmodeled  dynamies.  In  partieular,  it  was  shown  how  two 
properly  designed  filters,  one  in  the  eontrol  loop  and  one  in  the  parameter  estimation  bloek, 
yield  a  transient  response  due  to  adaptation  with  an  arbitrarily  small  maximum  norm.  Also, 
based  on  a  number  of  simulations,  it  was  verified  that  the  system  remained  stable  in  the 
presenee  of  external  disturbanees  and  unmodeled  dynamies. 

This  approaeh  was  extended  to  non-minimum  phase  systems  by  assuming  some 
stability  eonditions  using  a  non-proper  eompensator,  sueh  as  a  PID  eontroller.  In  the  latter 
ease,  the  effeetiveness  of  this  approaeh  was  demonstrated  by  applieation  to  a  system  with 
flexibilities,  having  a  non-minimum  phase  transfer  funetion. 

The  eontributions  in  this  dissertation  are  as  follows: 

1.  For  minimum  phase,  SISO  systems  global  stability  has  been  established 
together  with  eriteria  on  the  ehoiee  of  the  various  Alter  parameters.  Also, 
upper  bounds  of  the  transient  response  due  to  adaptation  have  been 
eomputed  and  shown  how  they  are  affeeted  by  the  filters’  bandwidths; 

2.  For  non-minimum  phase  systems,  we  have  shown  that  the  use  of  an 
augmenting  system  in  parallel  with  the  plant  makes  the  overall  system 
minimum  phase,  so  that  the  Adaptive  Controller  ean  be  applied.  The 
eriterion  of  ehoosing  the  augmenting  network  is  by  the  inverse  of  a  PID 
eontroller,  whieh  stabilizes  the  plant.  Although  the  PID  eontroller  has  a  non¬ 
proper  transfer  funetion  due  to  the  derivative  aetion,  its  inverse  is  proper 
thus  easily  implementable  in  the  presence  of  additive  disturbances.  The 
augmenting  network  is  chosen  so  that  it  has  no  effect  in  steady  state,  only 
during  the  transient  response  of  the  system; 

3.  Application  to  a  very  lightly  damped  flexible  system  with  two  degrees  of 
freedom  on  a  plane  shows  the  effectiveness  of  this  approach  in  an  actual 
implementation. 
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B,  FURTHER  WORK 

The  study  in  this  dissertation  presents  a  preliminary  investigation  on  an  Adaptive 
Control  system  with  the  intent  of  making  it  attraetive  to  engineering  problems.  There  are 
still  a  number  of  issues  that  need  to  be  addressed,  namely: 

1 .  Robustness  in  the  presenee  of  external  disturbanees  and  unmodeled 
dynamies.  Although  simulation  results  show  that  this  eontroller  still  keeps 
its  stability  properties  with  modeling  errors,  a  need  to  establish  how  mush 
uneertainties  ean  be  tolerated  while  still  guarantee  stability  has  to  be 
addressed.  In  partieular  the  elaim  of  arbitrarily  small  transient  errors  due  to 
adaptation  has  to  be  reformulated  to  take  modeling  errors  into  aeeount.  The 
expeeted  final  result  will  be  a  usual  tradeoff  between  performanee  and 
robustness; 

2.  The  extension  of  this  algorithm  to  non-minimum  phase  systems  seems  to  be 
very  important  in  view  of  the  faet  that  systems  with  lightly  damped 
flexibilities  are  very  hard  to  eontrol.  The  algorithm  presented  assumes  that  a 
stabilizing  non-proper  eontroller  (sueh  the  PID)  is  known  to  the  designer. 
The  issue  to  be  further  investigated  is  the  best  way  of  determining  this 
stabilizing  eontroller,  in  view  of  the  given  uneertainties  of  the  system  to  be 
eontrolled.  From  an  engineering  standpoint  this  would  eorrespond  to  the 
eombination  of  two  eontrollers:  a  standard  fixed  PID  eontrol  just  to  stabilize 
the  system  and  the  adaptive  eontroller  to  improve  its  dynamie  performanee. 
These  two  eontrollers  (fixed  and  adaptive)  interaet  with  eaeh  other,  but  their 
effeets  on  the  system  performanee  and  robustness  need  to  be  further 
investigated; 

3.  Finally  most  systems  of  interest  have  Multiple  Inputs  and  Multiple  Outputs 
(MIMO).  This  is  the  ease  of  the  flexible  arm  presented  in  the  last  seetion: 
extension  to  an  aetual  implementation  in  spaee  would  require  dynamies  in 
all  three  dimensions  with  aetuators  in  all  three  axes.  For  this  sort  of 
problems  there  are  two  possible  ways  of  extending  the  theory  and  design. 
One  way  is  to  take  the  whole  dynamies  of  the  system  into  aeeount  and 
extend  the  SISO  approaeh  to  MIMO.  Along  these  lines  a  eonsiderable 
amount  of  work  has  been  done  for  standard  Adaptive  Control  systems  and  it 
ean  be  extended  to  the  LI  approaeh.  However  in  some  eases  as  in  the 
flexible  arm  in  spaee,  the  MIMO  system  ean  be  modeled  as  a  set  of  SISO 
systems  with  eross  eoupling  dynamies.  In  this  ease  it  might  be  easier  to 
design  separate  independent  SISO  eontrollers  and  treat  the  eross  eouplings 
as  disturbanees.  In  this  way,  if  the  SISO  eontrollers  have  suffieient  degree 
of  robustness  the  overall  system  will  still  be  stable.  These  show  important 
tradeoffs  of  this  elass  of  design. 
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APPENDIX  A. 


Equation  (3.31)  can  be  written  as 


xxt) 


a 


-2  I - 


w,{t)\  +||X(0|| 


X^t) 


cr  w^{t)  +  X{t) 


-e{t) 


where  ju-jLi.{t)>0. 


(A.1) 


Define  the  positive  define  funetion 


v{e{t))  =  U{tye{t)  (a.2) 

then 

vm))  =  e\tym  (A.3) 


and 


vm))=-M- 


e\t)x{t) 


2N 


w^it)^+  X(t) 


-e(0  +  X(A-A,(0)- 


A,(0 


cr 


-2  I  - 


w,(0|  +||2'(<)|| 


■e(t).  (A.4) 


The  rightmost  term  is  nonzero  only  when 


4(0  =  ^  <  0  and  X.{t)e{t)  >  0 


(A.5) 


or 


L  4(0  =  ^-  ^  0  and  X,  (0e(0  <  0 

whieh  makes  it  always  non-positive. 

Substituting  for  the  error  term  e{t)  =  0^  {t)X{t)  +  w,,{t)  the  following  expression  is 
obtained 


V{0{t))  <  +s{t)ri{t)) 


-A 


^(0  +  ^?7(0 


(A.6) 
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where 


r 


s{t)  = 


e\t)x{t) 


< 


77(0  = 


+\\X{t)\\ 

^Eit) 


.s{t)  =  ■ 


e\t)X{t) 


V. 


^cr^|w^(0|^+||X(0f  ^a^\wE{t)f  +\\X{t)\ 


Now,  with  any  t,  T  >  0  and  integrating  Equation  (A. 6)  between  t  and  t+T 


(A.7) 


t+T 

0<vm+T))<v{e{t))-M\ 


t+T 

dT+—  J  \ri{T)'^  dr  (A.8) 


yields 


t+T 


dr  < 


2ju  4 


t+T 


j  dx. 


(A.9) 


Definition  of  r]{t)  in  Equation  (A.7)  yields  the  upper  bound  on  the  rightmost  term  as 


t+T  t+T 

j  dT<^  - 


0  if  w^(r)  =  0  for  all  r 


(A.10) 


We{t)\  j 

dx<  ^ 


otherwise. 


cr 


The  final  results  beeome 

a.  Equation  (3.33a),  eome  easily  sinee  the  initial  estimates  are  within  the 
bounds 

^  d{t)  <  ,  and 


4(0  4(0  >0 


(A.  11) 
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b.  From  the  definition  of  s{t)  in  Equation  (A.7) 


1 


(A.  12) 


where  we  ean  bound 


i, 


1  r 


d.T 


+ 


i: 


\Hr) 


1 

2  ^2 

dz 

) 


(A.13) 


Finally,  eombining  Equations  (A.9),  (A.12),  and  (A.13)  yields  Equation  (3.33b) 
with  bound  Equation  (3.34). 

e.  Sinee 


\\X{t)f\e{tf 


0{t) 


<A 


<A 


e{ty 


(a'|w^(0|'+||A(0f)  ^>^(01  +||A(0 


(A.  14) 


and  substituting  for  e{t)  =  6^  {t)X{t)  +  ^^(t)  yields 


e{t) 


yja^\w^{t)f  +\\X(t)\\ 


(  1^1 

=  s{t)  +  ri{t)=  s{t)  +  -ri{t)  +-?7(0 
V  1)1 


(A.15) 


This  expression  has  the  same  upper  bound  as  Equation  (A.13)  so  that  Equation  (3.33o) 
follows  easily  with  bound  again  of  Equation  (3.34) 


QED 
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APPENDIX  B. 


Proof:  From  Equation  (3.20),  the  expression  for  the  plant  output  y{t)  is 


E(0  = 


-J—  {kji (0  +  (^0  -  )u (0  +  h^(f)Xt)  +  k^ (0)  +  w{t). 
D{p) 


(B.l) 


Combining  Equations  (3.21),  (3.22),  and  (3.23),  and  rearranging  Equation  (3.20)  as 
follows  yields 


6’^^(0  =  D{p){y{t)  -  w{t))  +  (^0  - kju{t)  - k^uit).  (B.2) 


Erom  Equation  (3.1) 


-  >^(0) 


(B.3) 


Substituting  Equation  (B.3)  into  Equation  (B.2)  results  in 


=  D{p)y{t)  -  D{p)w{t)  +  (^0  -  k^ )u{t)  -  y{t)  +  (t). 

B(p)  B(p) 


(B.4) 


hen 


e^m  =  D{p) 


\  Ap)  ^ 

V  B{p)D{p)j 


y{t)  +  (^0  -  k^ )u{t)  -  D{p)w{t)  +  Wo  (t). 

B{p) 


(B.5) 


Substituting  Equation  (B.5)  into  Equation  (B.2)  yields 


( 


D{p) 
{^-C,{p)) 
D{p) 


A(p) 


V  B{p)D{p) 


yit) 


iK  Cp  {p))  w{t)  +  ^0  (0  +  w{t) 

D{p)  A{p) 


(B.6) 
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k^u{t)  =  Cp  {p)[e^  {t)(l){t)  +  v(0)  -  Cp  {p)D{p) 


1-- 


A{p) 


B{p)D{p) 


y{t) 


A(p) 


-Cp  (p)(kg  -  )u{t)  +  Cp  {p)D{p)w{t)  -  Cp  {p)  ——  Wo {t). 

B{p) 


(B.7) 


Equation  (B.6)  can  be  adjusted  as 
y{t)  =  [\-Cp{p)) 


\  Ap)  ^ 

V  B{p)D{p)j 


E(0  +  {0^ +  v{t)) 


-w{t)  +  Cp  {p)w{t)  +  [\-Cp{p)) 


D{p) 

Ap) 

yB{p)D{p)j 


D{p) 


wA)  +  w{t) 


then 


f 

y(t)  =  {\-Cp(p))  1 

V 


Ap) 

B{p)D{p)^ 


y(t)  + 


i^-Cpip)) 

Dip) 


{rKAt)) 


+ 


Dip) 


+  {^-Cpip)) 


Ajp) 

Bip)Dip)^ 


wA)  +  Cpip)wit) 


(B.8) 


and  finally 


(1  +  rCp  ip))  At)  =  -Cp  ip)Dip) 


1-- 


Aip) 


Bip)Dip)^ 


-Cpip) 


Ajp) 

Bip) 


Woit)-Dip)wit) 


yit)  +  Cp  ip)[A  it)  At)  +  v(0) 

(B.9) 


Then  k^pit)  is  isolated; 


KAt)  = 


-Cpip)Dip) 


f 


Aip) 


{\  +  yCpit)) 

Cpip)  (Ajp) 
{\  +  yCpit))yBip) 


V  Bip)Dip) 

f 


wA)-Dip)wit) 


(B.IO) 


Substituting  kjiit)  from  Equation  (B.IO)  into  Equation  (B.8),  results  in  the  proof  of 
Equations  (3.41)  and  (3.42). 
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where 


y{t)  =  H„{p)y{t)  +  [p)Cp{p){e^  {t)(l){t)  +  v(0)  +  H^[p)  w^{t) 

kji{t)  =  {p) y(t)  +  {p)Cp(p)[0^ (t)<p(t)  +  v(0)  +  //„„  (p)  Wo(t) 


(3.41) 


and  / 


Hjp)  =  {i-CAp)) 


B(p)D(p) 


\+rCp{p) 


(i^rC,{p))Dip) 

H  (n\-  I  I 

{\  +  yC^{p)){B{p)D{p)  P{p)D{p)  j 


^UO(P) 

^uy(p) 

Hu«{p) 


1 

\  +  YCp{p) 

Cf[p)D{p)(  A{p) 

\  +  YCp[p)  yB{p)D{p)  ^ 

Cp{p)D{p)(  A(p)  ^  kjp) 

\  +  YCf{p)  yB{p)D{p)  P{p)D{p)  ^ 


kjp) 

P{p)D{p) 


k  k 

=  —  - 1  in  the  interval  0  <  /  <  Ymax  -  “  1  • 

k^  k^ 


(3.42) 


QED 
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APPENDIX  C. 


Let 

a.  C^{p)  be  the  transfer  function  of  an  exponentially  stable,  causal  Linear  Time 
Invariant  (LTI)  system, 

b.  0(t)  be  a  bounded  vector 

c.  (p(t)  e  ,e{t)  e  i?  be  the  input  and  output  of  the  following  Linear  Time  Varying 
operator 


~e{t) = c,  {p)[e^  im))  -  (OQ  {p)m) 


(C.l) 


then 


t  t 

e{t)  =  I {T)(l){T)dT  {t)(l){T)dT  (C.2) 

0  0 


Figure  C.L  Region  of  Integration  for  Equation  (C.2). 
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Since 


e{t)  =  e{Qi)  +  \^{X)dX 
0 


Substituting  Equation  (C.3)  into  Equation  (C.2)  yields 


(C.3) 


t  t 

e{t)  =  -T)6{X)(l){T)dMT  (C.4) 

0  T 


It  is  easy  to  see  from  Eigure  C.  1  that  the  region  of  integration  can  be  expressed  as 


0<X<t 

0<T<X 


(C.5) 


which  yields 


t  .  i 

e{t)  =  -J  0^  (/l)|  {t  -  z)(j){T^dTdX 
0  0 


A  simple  change  of  integration  variables  yields  the  result. 


(C.6) 


QED 
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